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Abstract 



The central theme of this paper is the variational analysis of homeo- 
morphisms h: X Y between two given domains X, Y C M". We look for 
the extremal mappings in the Sobolev space >^^'"'(X, Y) which minimize the 
energy integral 



Because of the natural connections with quasiconformal mappings this n- 
harmonic alternative to the classical Dirichlet integral (for planar domains) 
has drawn the attention of researchers in Geometric Function Theory. Ex- 
plicit analysis is made here for a pair of concentric spherical annuli where 
many unexpected phenomena about minimal n-harmonic mappings are ob- 
served. The underlying integration of nonlinear differential forms, called free 
Lagrangians, becomes truly a work of art. 
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Preface 



The future developments in modern geometric analysis and its governing 
partial differential equations (PDEs) will continue to rely on physical and 
geometric intuition. In recent years, this trend has become more pronounced 
and has led to increasing efforts of pure and applied mathematicians, en- 
gineers and other scientists, to share the ideas and problems of compelling 
interest. The present paper takes on concrete questions about energy min- 
imal deformations of annuli in R" . We adopted the interpretations and 
central ideas of nonlinear elasticity where the applied aspects of our results 
originated. A novelty of our approach is that we allow the mappings to slip 
freely along the boundaries of the domains. It is precisely in this setting 
that one faces a real challenge in establishing the existence, uniqueness and 
invertibility properties of the extremal mappings. The underlying concept 
of Free Lagrangians, is the core of the matter. 

Our approach is purely mathematical though the questions are inti- 
mately derived from Nonlinear Elasticity. Both the theoretical and prac- 
tical aspects of this work culminate in actual construction of the mappings 
with smallest conformal energy. Special efforts have been devoted to some- 
what subtle computational details to present them as simply and clearly as 
possible. 

We believe the final conclusions shed considerable new light on the Cal- 
culus of Variations, especially for deformations that are free on the boundary. 

We also feel that some new facts discovered here have the potential for ap- 
plications in Geometric Function Theory as well as for better understanding 
the mathematical models of Nonlinear Elasticity. 
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CHAPTER 1 



Introduction and Overview 



1.1. Basic notation 



Let us take a moment to recall a very much needed notation from the cal- 
culus of vector fields and matrix fields. Consider a mapping h : X '-^ Y be- 
tween domains X C M" and Y C M"", h = {h^, h?,..., N^), where h^,...,N^ 
arc scalar functions in the Sobolev space W^f(%.). The differential Dh{x), 
defined at almost every x G X, represents a linear transformation of M" into 
M™, Dh{x) : MP M™. With the standard choice of the coordinates in 
and we have a matrix field, again denoted by Dh, 



(1.1) 



Dh 



hi h 



X2 



hi 



hi 



Hereafter, we abbreviate the notation of distributional partial derivative 
^ to Fx^, i = l,...,n for F G ^j^^,(X). The differential matrix, also 
called Jacobian matrix or deformation gradient, acts on a vector field V = 
{V\V^,..., y") G ifi^^(X,M") by the rule 



Here V stands for the gradient operator acting on real-valued functions in 
#^^^^(X). More generally, consider an arbitrary matrix field 





' h^ 


■ hi 








' (V/ii,F) " 


[Dh]V = 
















"■xi "'X2 


• K 




yn 




_ (V/i™,y) 



(1.2) 



M 



Ml{x) 



Mf (x) 



M™(x) 



and denote its row- vector fields by r^, . . . ,r™ G ^j^^(X, M"). Similarly, the 
column-vector fields will be denoted by ci, . . . , c„ G ^^^(K,W^). The di- 
vergence operator acting on a vector field r = (ri, . . . ,r„) G ^^^(K,M.'^) is 
a Schwartz distribution defined by 
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Then the divergence of a matrix field M e 
valued in W^, 

divr^ 

DivM 



') is a distribution, 



divr" 



In particular, 



where A = ^ + • • • + is the usual Laplacian. The matrix fields M G 

^^p(X, M"*^") which satisfy the equation DivM = will be called divergence 
free, meaning that 

/ (M, Dri) = for every test mapping r/ G ^o°°(X, R"*). 

Hereafter {A,B) = Tr{A*B) is the inner product of matrices. We will be 
typically working with the Hilbert-Schmidt norm of a matrix 

m n 

\\M\\' = {M,M)=J2J2\^if- 

j=l i=l 



1.2. Mathematical Model of Hyperelasticity 

Geometric Function Theory (GFT) is currently a field of enormous ac- 
tivity where the language and general framework of Nonlinear Elasticity is 
very helpful. As this interplay develops, the n-harmonic deformations be- 
come well acknowledged as a possible generalization of mappings of finite 
distortion. We have also found a place for n-harmonic deformations in the 
theory of nonlinear hyperelasticity. J. Ball's fundamental paper [5] accounts 
for the principles of this theory and sets up mathematical models. Histor- 
ically, the relation between hyperelasticity and quasiconformal theory, has 
not been clearly manifested, but it is indeed very basic and fruitful. 

One can roughly describe the hyperelasticity as a study of weakly dif- 
fercntiablc homcomorphisms h: X Y between domains in M" (or n- 
manifolds) that minimize a given energy integral. 



(1.3) <4 = [ B{x,h,Dh)dx <oo D/i:X^M'*^", 

The condition on the injectivity of h is imposed in order to avoid interpen- 
etration of matter. The Jacobian matrix Dh{x) G M"^" , defined at almost 
every point a: G X , is referred to as the deformation gradient. In this model 
the so-called stored energy function E: X x Y x M"^" — )• M is given; it 
characterizes mechanical properties of the elastic material in X and Y . 
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Motivated by GFT we will be essentially concerned with the n-harmonic 
energy, also called conformal energy 

(1.4) <^h= [ \\Dh{x) irdar. 



Another energy integral of interest in GFT is 

(1.5) 



We devote Chapter 8 to this latter integral. 

1.3. Variational Integrals in GFT 

In another direction, we recall Geometric Function Theory in R" and 
its governing variational integrals. Let as begin with a conformal mapping 
h:X^Y. Thus at every a; G X we have the relation between the norm of 
the Jacobian matrix and its determinant || Dh{x) \\ " = J[x, h). This can 
be expressed in the form of a nonlinear Cauchy-Riemann system of PDEs; 



(1.6) D*h- Dh = J{x,h)n I 



It is evident that the n -harmonic energy of /i : X Y depends only on 
the deformed configuration. Indeed, we have 



n2|Y| 



(1.7) ^h= [ II Dh{x) r dx = /" J{x, h) dx = 

For other homeomorphisms g : X Y, in the Sobolev space #'^''^(X, Y) , 
we only have a lower bound, due to Hadamard's inequality for determinants: 
r r 
Sg= I II L»5(a;) IP dx > nt / J{x,g)dx = n'^\Y\ 
Jx Jx 

Thus conformal deformations h : X Y are none other than those hav- 
ing the n-harmonic energy equal to n^jY], the smallest possible. It is for 
this reason that conformal mappings are frequently characterized as absolute 
minimizers of the n-harmonic integral. However, it is rare in higher dimen- 
sions that two topologically equivalent domains are conformally equivalent, 
because of Liouville's rigidity theorem. Even in the plane, multiply con- 
nected domains like annuli are of various conformal type. From this point 
of concerns Quasiconformal Theory [3, 29] offers significantly larger class of 
mappings. 

Definition 1.1. A homeomorphism /i : X — > M" of Sobolev space 
>^^^^(X, R") is said to have finite outer distortion if 

II Dh{x) ir ^ n"^ K{x) J{x,h) 

for some measurable function 1 ^ K{x) < oo. The smallest such K{x) is 
called the outer distortion, denoted by [x, h) . Then h is K - quasiconformal 
if Kq {x, h) ^ K for some constant K. 
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A concept somewhat dual to outer distortion is the inner distortion. For 
this, consider the cofactor matrix D'^h (it represents infinitesimal deforma- 
tions of (n— l)-dimensional entities) defined for invertible Dh via Cramer's 
rule 

D^h- D*h = J{x,h)I 

For a map h G >^^^^(X, M'^), not necessarily a homeomorphism, but with 
nonnegative Jacobian, we introduce the inner distortion function 

(1.8) KJx,h) = \ ""/^ [J(-,h)]-^ 11 J[x,n)>u 

\l if J{x,h)=0 

Remark 1.2. Any map of finite outer distortion has finite inner distor- 
tion and ]Kj-(x) < K^{x)"'~^ , but not vice versa. 

It is again interesting to find the place for such mappings in contin- 
uum mechanics. The latter deals with the positive definite matrix C{x) = 
D*h{x) Dh[x) as the right Cauchy-Green deformation tensor. While on the 
other hand, there is a fundamental interplay between mappings of finite 
distortion and the Beltrami equation 

(1.9) D*h{x)Dh{x) = J{x,h)n G{x) , det G(a;) = 1 

Thus G = G{x), called the distortion tensor of h, is none other than the 
Cauchy-Green tensor renormalized so as to have determinant identically 
equal to one. The symmetric positive definite matrix function G = G(a;) = 
[Gij (x)] G M"^" can be viewed as a Ricmann metric tensor on X. In this way 
h becomes conformal with respect to this, usually only measurable, metric 
structure on X. Thus G(x) is uniformly elliptic in case of iC-quasiconformal 
mappings. 

It is in this Riemannian manifold framework that variational interpreta- 
tions of quasiconformal mappings really crystalize. For example, the solu- 
tions to the Beltrami equation (1.9) are none other then the absolute mini- 
mizers of their own energy integrals. Indeed, a homeomorphism /i : X — t- Y 
of Sobolev class >^^''*(X, Y) solves the Beltrami equation (1.9) if and only if 

(1.10) (4= / B{x,Dh)dx = nt [ J{x,h)dx = nt|Y| 

Jx Jx 

where the integrand is defined on XxR"^" by the rule 

n 

E(x,0= {Tr[^G-\x)e]y , 

As in the conformal case, for all other homeomorphisms g: X Y, in 
the Sobolev space /^^'"(X, Y), we have the lower bound S'g ^ n2[Y]. The 
most appealing conclusion is a connection between the n -harmonic energy 
of : X — )• Y and the inner distortion function of the inverse mapping 
f = h-^ : Y ^ X. 
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Proposition 1.3. (Transition to the inverse map) 
Let f G /^^^""^(Y, X) be a homeomorphism of finite outer distortion between 
bounded domains, with IC^(y, /) G ^^(Y). Then the inverse map h = f"^ : 
X Y belongs to the Sobolev class >^^'"(X, Y) and we have 

(1.11) nt / K^iy,f)dy = [ WDhixWda: 

This identity gains in significance if we realize that the polyconvex vari- 
ational integrand in the left hand side turns into a convex one, a rarity 
that one can exploit when studying quasiconformal mappings of smallest 

-mean distortion. From yet another perspective, it is worth mentioning 
the classical Teichmiiller theory which is concerned, broadly speaking, with 
extremal mappings between Riemann surfaces. The extremal Teichmiiller 
mappings are exactly the ones whose distortion function has the smallest 
possible ^°^-norm. The existence and uniqueness of such an extremal qua- 
siconformal map within a given homotopy class of quasiconformal mappings 
is the heart of Teichmiiller 's theory. Now, in view of the identity (1.11), 
minimizing the ^^-norm of the inner distortion function offers a study of 
n-harmonic mappings. Is there any better motivation? 

1.4. Conformal Energy 

For X and Y open regions in R", we shall consider mappings 

(1.12) h= (/i\...,/i") : X — ^ Y 

in the Sobolev class W^''^{X,Y). Thus the Jacobian matrix of h and its 
determinant are well defined at almost every point x & X. We recall the 
notation 

G=^"(X,M"^"), J{x,h) = detDh 

Here, as usual, M"^" is supplied with the inner product and the Hilbert- 
Schmidt norm: 

n 

(1.14) {A, B) =Tr {A*B) = ^ A)B'j \\A\\ = {A, A)h 

ij=l 

At the initial stage of our undertaking the n-harmonic integral will be sub- 
jected to the orientation preserving homeomorphisms h : X Y , so 
that 

r r 

(1.15) <^h= \\Dh{x)\\''dx ^ nt / J{x,h)dx = ni|Y| 

Jx Jx 

In dimension n ^ 3, it may well be that no homeomorphism h : X Y of 
finite n-harmonic energy exists, as the following result [31] shows. 



(1.13) Dh 



dK_ 
dx 



3. 
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Theorem 1.4. Let X C M"^ he a hall with a k- dimensional closed disk 
removed, and let Y C he a hall with a {k + \)- dimensional closed disk 

removed, \ ^ k < n — \. Then every homeomorphism h : X Y has 
infinite n-harmonic energy. 

Note that both X and Y arc topological annuli; that is, homeomorphic 
images of a spherical annulus A = {x: r < \x\ < R}. Let us view the 
disks removed from the balls as cracks. It can be easily shown, by means 
of an example, that mappings of finite conformal energy may outstretch an 
(n— l)-dimensional crack into an n-dimensional hole. However, Theorem 1.4 
ensures us that, in principle, mappings of finite energy cannot increase the 
dimension of lower dimensional cracks, a fact highly nontrivial to observe 
and prove [31]. From now on we assume without explicit mention that the 

domains X and Y admit at least one homeomorphism h : X Y in the 
Sobolev space ^i'"(X,Y) 

1.5. Weak limits of homeomorphisms 

But the true challenge is to find a deformation h° : X Y with the 
smallest possible energy. In general, when passing to the weak limit of the 
minimizing sequence of homeomorphisms, the injectivity of the extremal 
map will be lost. Nevertheless, from the point of view of the elasticity the- 
ory [1, 5, 10], such limits are still legitimate deformations to consider. For, 
if this is the case, they create no new cracks or holes in Y. Let ,^(X, Y) 
denote the class of weak limits of homeomorphisms h : X Y in the 
Sobolev space W^''^(K,Y). We refer to such limits as permissible deforma- 
tions.^ Differentiability and geometric features of permissible mappings are 
not as clear as one may have expected. In Theorems 1.5, 1.6 and 1.7 we 
assume that X and Y are bounded domains of the same topological type, like 
spherical annuli, having at least two though finitely many boundary compo- 
nents. We consider a sequence hj : X Y of homeomorphisms converging 
weakly in W^''^{'K, Y) to a mapping /i: X — >■ Y. In general, homeomorphisms 
hj : X Y do not extend as continuous maps to the closure of X , but 
the distance functions x — )• dist(/ij(a;), OY) do extend. This is also true for 
the limit mapping h. The precise result to which we are referring is the 
following: 

Theorem 1.5. [30, Theorem 1.1] For the ahove-mentioned pair of do- 
mains X and Y, there exists a nonnegative continuous function rj = r]{x) 

defined on X such that 

(1.16) dist(/ij(x),5Y) ^ r]{x) \ Dhj \\ ifn(x)) rf = on dX. 



^Homeomorphisms converging weakly in #'^'"(X, Y) also converge c-uniformly, so 
their limits are still continuous, taking X into Y . 
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dist(/i(x), 5Y) ^ r]{x)sup \\Dhj || jfn(x)- 



The weak limit h actually covers the target domain, but this may fail if 
X and Y have only one boundary component. 

Theorem 1.6. [30, Theorem 1.4] The mapping h is continuous and 
Y C /i(X) C Y. Furthermore, there exists a measurable mapping / : Y — )■ X, 
such that 



everywhere on Y. This right inverse mapping has bounded variation, 



As noted in [30, Remark 9.1] the weak limit h is monotone in the sense 
of C.B. Morrey [37]. 

Theorem 1.7. The mapping h is monotone, meaning that for every 
continuum IfC C Y its preimage C X is also a continuum; that is 

compact and connected. 

The proof of this theorem is presented in Section 3.3. 



The first nontrivial case is that of doubly connected domains. Thus we 
consider mappings /i : A — t- A* between concentric spherical annuli in M". 

A = A(r, i?) = {x G M" ; r<\x\<R}, ^ r < R < oo 
A* = A(r. , i?J = {y G R" ; r^ < \y\ < R,} , ^ < i?, < oo 
Such domains are of different conformal type unless the ratio of the two radii 
is the same for both annuli. As for the domains of higher connectivity in 
dimension n = 2 , the conformal type of a domain of connectivity ^ > 2 is 
determined by M — 6 parameters, called Riemann moduli of the domain.^ 
This means that two ^-connected domains are conformally equivalent if and 
only if they agree in all 3i — 6 moduli. But we shall have considerably 
more freedom in deforming X onto Y, simply by means of mappings of 
finite energy. An obvious question to ask is whether minimization of the 
n-harmonic integral is possible within homeomorphisms between domains 
of different conformal type. 

Concerning uniqueness, we note that the energy S'^ is invariant under 
conformal change of the variable a:; G A. Such a change of variable is realized 
by conformal automorphism of the form 



o / = id : Y ^ Y, 



/ 



BV(Y) 




1.6. Annuli 



(1.18) 




in this context the mappings are orientation preserving. 
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where k = 0,1 and T is an orthogonal matrix. 

1.7. Hammering a part of an annulus into a circle, n = 2 

Let us caution the reader that a minimizer h° : A A*, among all 
permissible deformations, does not necessarily satisfy the Laplace equation. 
A loss of harmonicity occurs exactly at the points where h° fails to be 
injective. This is the case when the target annulus A* is too thin as compared 
with A; precisely, if 

R \ f R t\ 

(1.19) — - < - ( 1- t: I ! -annuli below the Nitsche hound 

r* 2 \ r RJ 

By way of illustration, consider the so-called critical Nitsche map 

(1.20) ^{z) = 1 + , < |z| < GO 

This harmonic mapping takes an annulus A(l,/?) univalently onto A* = 
A(l, R^,) , where R^ = \{^R + We have equality at (1.19), and N is the 
energy minimizer. Note the symmetry H (=) = '^{z). Thus the same Nitsche 
map takes reflected annulus A(i?~^,l) univalently onto A*. Let us paste 
these two annuli along their common boundary 

A = A(r,i?) = A(r,l] UA[l,i?) , r = - 

R 



Now the same harmonic map ^ : A[l,i?,) is a double cover. Its 

Jacobian determinant vanishes along the unit circle, the branch set of 
Therefore, K is not permissible (it is not a weak /^^'^-limit of homeomor- 

phisms). An extension of : A(l,i?) A* inside the unit disk to a 
permissible mapping of A(r, R) onto A* can be nicely facilitated by squeez- 
ing A(r, 1) onto the unit circle. This procedure will hereafter be referred 
to as hammering the inner portion of the domain annulus onto the inner 
boundary of the target. Precisely, the map we are referring to takes the 
form 



h%z) = H{\z\)^M. 



||| < |z| ^ 1 , hammering part 

^ ^ (z + I) 1 ^ |2;| ^ i? , harmonic part 

It is true, though somewhat less obvious, that: h° is a /^^'^-limit of 
homeomorphisms h: A A* and its energy is smaller than that of any 
homeomorphism from A onto A*, see Figure 1. 

We shall actually prove the following theorem. 

Theorem 1.8. Let A = A{r,R) and A* = A(r^,R^) be planar annuli, 
< r < R < oo and < r* < i?* < oo. We have: 
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R* 
1 




r 1 R t 



Figure 1. Hammering the inner ring into the unit circle. 



Case 1. (Within the Nitsche hound) If 



R* ^ 
~ ^ 2 

then the harmonic homeomorphism 



R r 
V'^R 



attains the smallest energy among all homeomorphisms h: A A*, and as 

such is unique up to a conformal automorphism of A. 
Case 2. (Below the Nitsche bound) If 



R^ \(R r_ 
r, ^ 2[r^ R 



A* is not 



a 

'a^R 



-critial Nitsche configuration. 



then the infimum energy among all homeomorphisms h: A 
attained. Let a radius r < a < R be determined by the equation 

Rif 1 / R 
rT ~ 2 ^ 

Then the following mapping 

{r^A r <\z\^a 

^(f + f) ^^\z\<R 

is a W^''^-limit of homeomorphisms hj : A A*, and its energy is smaller 
than that of any homeomorphism h: A A*. 

The proof of this theorem was first given in [2]; here in Section 10 we 
present another one, based on free Lagrangians. 

Let us emphasize that this result does not rule out the existence of 
univalent harmonic mappings from A onto A* , simply because harmonic 
homeomorphisms need not be the ones that minimize the energy. Nonex- 
istence of harmonic homeomorphisms between such annuli A and A* was 
conjectured by J. C. C. Nitsche [41]. After several partial results were ob- 
tained in [17, 34, 35, 50], the conjecture was proved in [25, 26]. The 
connection between the Nitsche conjecture and minimal surfaces is further 
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explored in [28]. Hammering also occurs in free boundary problems for min- 
imal graphs, where it is called edge-creeping [9, 21, 46]. Similar hammering 
phenomena will be observed in higher dimensions as well. 

1.8. Principal n-harmonics 

In studying the extremal deformations between spherical annuli it is nat- 
ural to look for the radially symmetric solutions of the n-harmonic equation 

(1.21) div ( II m ip-^m) = , h{x) = H{\x\)-^^ 

There is a nice reduction of this problem to the first order (nonlinear) dif- 
ferential equation for the strain function H = H{t), called the characteristic 
equation for H, 



(1.22) CH 



„2 



n-1 



n — 2 
2 



(^H"^ - t^ijsj = const. 



Although this equation provides a very convenient tool for studying proper- 
ties of radial n-harmonics, a little caution is needed because ^-solutions to 

(1.22) may fail to satisfy the original equation (1.21). We shall distinguish 
four so-called principal n-harmonics , ^° , and K_. The first two are 
conformal mappings 

(1.23) Ko(ar) = x and K°(x) = ^ 

The other two are more involved, but still radially symmetric 

(1.24) K+ = i/+(|a;|)^ and H_ = //_(|x|)^ 

They cannot be described in any elementary way, except for the case n = 2. 

(1.25) Ho{z) = z,H%z) = ^ 

(1.26) K+(z) = ^ + =) and K_(^) = 1(^-1 

The principal n-harmonics and N_ are determined by solving the follow- 
ing Cauchy problems for their strain functions 

(1.27) CH+ = 1 H+{1) = 1 

(1.28) CH^ = -1 iJ_(l) = 

It has to be emphasized that all radial solutions to the n-Laplace equation 
take the form 

(1.29) /i(x) = AH(fex), AgR,/c>0 

where is one of the four principal solutions. We then sec that every radial 
solution (originally defined in an annulus) extends n-harmonically to the 
entire punctured space = M'* \ {0}, and is evidently "^""-smooth. Their 



1.9. ELASTICITY OF STRETCHING 11 

connections with (1-26) motivate our calhng Ho , , and the Nitsche 
maps in M". 

1.9. Elasticity of stretching 

There is an important entity associated with the radial mappings; namely 
the elasticity of stretching 

(1.30) ri„ (t) = , provided + ^ 

The elasticity function of a power stretching H{t) = t° is a constant, equal 
to a. Two particular cases ry = 1 and r] = —1 correspond to the conformal 
maps Ho and In both cases CH = 0. There are exactly three types of 
radial n-harmonics: 

(a) Conformxilly Expanding. These are the mappings with |ry^| > 1, ev- 
erywhere. Equivalently, CH = const > 0. Geometrically it means 
that h exhibits greater change in the radial direction than in spher- 
ical directions. 

(b) Conformally Contracting. These arc the mappings with |7y^| < 1, 
everywhere. Equivalently, CH = const < 0. 

(c) Conformally Balanced. These are conformal mappings with |?7^| = 
1, or, equivalently CH = 0. 

For every permissible radial map h{x) = H(\x\^-^^ : A — )■ A* (weak W^'""- 
limit of homeomorphisms) the elasticity function does not change sign. In- 
deed, this follows from the identity J{x,h) = H{\x\) ■ ^^^^^^ ■ The 

following ^^-mean (with respect to the conformal density on A) is equal 
to the ratio of the moduli of the annuli, see (1.40) for the definition of the 
modulus, 

, , /" I / M , / X Mod A* , dx t 

(1.31) -4 |r?„(x)| dfi{x) = , d^ 



Mod A 

In this way we are led to three types of pairs of the annuli A and A*. 

(1) Conformally expanding; it pertains to a pair of annuli such that 

Mod A* > Mod A 

(2) Conformally contracting; it pertains to a pair of annuli such that 

Mod A* < Mod A 

(3) Conformally equivalent; these are the annuli having the same mod- 
ulus 

Mod A* =ModA 

These three cases will be treated by using somewhat different estimates. 



"^The integral mean notation -ff{x)dij,{x) stands for the ratio J^fd^/ /j-d/tt. 
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1.10. Conformally expanding pair 

We shall show that the principal solution = generates 

all minimizers of the conformal energy. It is rather easy to show that for 
a given expanding pair A, A* there exist unique A; > and A > such 
that the n-harmonic map h°{x) = AH_(A;x) takes A homeomorphically onto 
A*. However, the answer to the question whether this map minimizes the 
n-harmonic energy among all homeomorphisms is not obvious. When n = 2 
or n = 3, the answer is "yes". 

Theorem 1.9. Let Mod A* > Mod A. Then for n = 2,3, then-harmonic 
radial map h° = X'ii-{kx) assumes the minimum conformal energy within 
all homeomorphisms. Such a minimizer is unique up to a conformal auto- 
morphism of A. 

Surprisingly, for n ^ 4 the answer will depend on how wide is the target 
annulus A*, relatively to A. 

Theorem 1.10. For dimensions n ^ 4, there exists a function J\f^ = 
J\f^{t), t <J\f^t) < oo for t>0, see Figure 2, such that: if 

(1.32) Mod A ^ Mod A* ^ Art(Mod A) -upper Nitsche bound for n ^ 4, 

then the map /; ° : A — )• A* is a unique ( up to an automorphism of A) 
minimizer of the conformal energy among all homeomorphisms. 

Theorem 1.11. In dimensions n ^ 4, there are annuli A and A* such 
that no radial stretching from A onto A* minimizes the conformal energy. 

1.11. Conformally contracting pair 

In this case we obtain the minimizers from the principal solution = 
i?_l_(|a;|) As before, we observe that the mappings 

(1.33) h°{x) = \H+{kx), k>Q, A>0 

are radial n-harmonics. Recall that A* is conformally thinner than A. But it 
is not enough. In contrast to the previous case, such n-harmonic mappings 
take the annulus A homeomorphically onto A* only when A* is not too thin. 
The precise necessary condition reads as 

(1.34) the lower Nitsche bound; A/t(ModA) ^ Mod A* ^ Mod A. 
Numerically, the lower Nitsche function J\f\ is given by 

(1.35) 0<^{t)=uJn-i^ogH+(exp—!—]<t for < t < oo. 

Theorem 1.12. Under the condition at (1.34) there exist unique k > 
and A > such that h°{x) = XH^{kx) takes A homeomorphically onto A*. 
This map is a unique (up to a conformal automorphism of A) minimizer of 
the conformal energy among all homeomorphisms of A onto A* . 



1.12. THE CONFORMAL CASE Mod A = Mod A* 
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Rather unexpectedly, the injectivity of the weak Umit of a minimizing 
sequence of homeomorphisms hj : A A* will fail once the lower bound 
in (1.34) is violated. Let us look more closely at the critical configuration 
of annuli; that is, 

(1.36) A4(ModA')=M„dA-, { I^r^fl. 2 

Thus the critical Nitsche map H+(a;) = H^{\x\) j^, defined for 1 ^ |x| ^ i? 
takes A' homeomorphically onto A* . The Jacobian determinant of K vanishes 
for |x| = 1. Now let us built a pair of annuli below the lower bound at (1.34), 
simply by pasting an additional spherical ring to A' along the unit sphere. 
Thus, we consider a slightly fatter annulus 

A = A{r,R)= A(r, 1] U A[l, R) 

As in dimension n = 2, we have 

Theorem 1.13. The following deformation 

. 0/ X def f ifi r < |x| ^ 1 , hammering part 

(1.37) (2^) = wf / \ II 7-. 7 

[ H+(x) 1 ^ |x| ^ it , n-harmomc part 

is a W^'"^ -limit of homeomorphisms h: A A* and its energy is smaller 
than that of any homeomorphism from A onto A*. Among all such mappings 
h° is unique up to conformal automorphisms of A. 

1.12. The Conformal Case ModA = ModA* 

Obviously the minimizers are conformal mappings. Precisely, they take 
the form 

±1 



(1.38) h°{x) = ^^^^ ( ^ 1 Tx 

\x\ \ \x\ I 

where T is an orthogonal transformation. This case receives additional con- 
sideration in Section 15. 
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1.13. The energy function 

So far we considered the domains X and Y equipped with the Euchdean 
metric. However, one may ask analogous questions for different metrics 
on X and Y. Wc do not enter into a general framework here, but instead 
illustrate this possibility by introducing a conformal density on the target 
space, which we continue to assume to be a spherical annulus 

dy 



(1.39) Y = A* = A(r* , i?*) , equipped with the measure d/x(y) 

The modulus of A* is none other than its "conformal" volume 

(1.40) ModA*=/ dii{y)= f i^=a;,_ilog:^ 

JK* J A* \y\ f* 

More generally, \i h : X A* is conformal then the modulus of X is 
defined by pulling back the measure d/i to X. 

J{x,h)Ax _n f \\Dh{x)\\'^dx 



(1.41) ModX= /%^=--^ / 



\h{x) 



Obviously, this definition is free from the choice of the conformal map h . 
For n = 2 the modulus of a doubly connected domain is the only conformal 
invariant; that is, the Riemann moduli space is one dimensional. The situ- 
ation is much more rigid for n ^ 3. If now h : X A* is any permissible 
deformation, then 

Dh(x) W^dx n f J(x,h)dx , 
V ^ II ^ / V > ; _ n2ModA* 



/, io\ a?- def f \\Dh{x)\\"-dx n f 

1.42 = / |w Vn / 



\h{x)\n 



Naturally, this integral tells us about how much h differs from a conformal 
mapping in an average sense. On the other hand, Quasiconformal The- 
ory deals with point-wise distortions. Among them is the outer distortion 
function 



(1.43) /Co(x,/i) = <{ n'^JixA) ' 

[ 1 , otherwise 

Let us push forward /Cq to the target space via the mapping h itself, so as 
to obtain a function 

(1.44) }Ch{y) = IC^{x,h), where x = h-'iy) 
We note, without proof, that 

(1.45) ^h = n^ f ^h{y)My) 

Jy 

This discussion leads us to the minimization of the =Sf ^(Y, d/i)-norm of the 
outer distortion. One might suspect that the minimum will be attained 



ij^mii^ ifj(x,/i)>o 



1.14. FREE LAGRANGIANS 
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when the distortion function is constant. Indeed, when the domain is also 
an annulus, say X = A = A(r, R), then the power stretching 
(1-46) ^^^^^ 

h°'(x) = A|a;|°~^x , where a = — — and A = r*r~" = R*R~°' 

^ ' ' ' ModA 

has constant outer distortion /Cq(x, h"^) = a~^n~^ (a^ + n — 1)?. Neverthe- 
less, it takes some effort to show that in dimensions n = 2, 3, /i" is actually 
a minimizer of among all homeomorphisms. 

Theorem 1.14. Let A and A* he spherical annuli in M", n = 2, 3. Then 
for every homeomorphism h : A A* we have 

(1.47) I i^^^^(n-l + a2)^ModA, where a = 

Ja l^h ModA 

Equality holds for the power stretching h{x) = r* r~'^\x\°'~^x, uniquely up to 
conformal automorphisms of A. 

It is somewhat surprising that in dimensions n ^ 4, this feature of no 
longer holds when the target annulus is conformally too fat. We have the 
following result. 

Theorem 1.15. For each n ^ A, there exists a„ > 1 such that (1-4V 
holds whenever 

n def ModA* 

«=m;^<"" 

The power stretching h{x) = r*r~° is the only minimizer of 
modulo conformal automorphism of A. 

Examples will be given to show that the extremals are no longer power 
stretchings if 



(1.49) I" 



ModA V n-3 

In other words, the upper bound for a in (1.48) lies in the interval 1 < 

an < \J~^^- See Chapter 12 for more precise estimates of an- Moreover, if 
ModA* is too large relative to ModA, then the extremals cannot be found 
even within general radial mappings, see Chapter 14. 

1.14. Free Lagrangians 

In 1977 a novel approach towards minimization of polyconvex energy 
functionals was developed and published by J. Ball [5]. The underlying 
idea was to view the integrand as convex function of null Lagrangians. The 
term null Lagrangian pertains to a nonlinear differential expression whose 
integral over any open region depends only on the boundary values of the 
map, like integrals of an exact differential form. The interested reader is 
referred to [6, 15, 22]. But we are concerned with mappings /i : X ^ Y 
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that are free on the boundary. The only condition we impose on h is that it 
is a weak /^^'"-hmit of homeomorphisms from X onto Y. There still exist 
some nonlinear differential forms, associated with a given pair of domains X 
and Y, whose integral means over X remain the same within a given class 
of deformations h: X Y, regardless of their boundary values. These 
are rather special null Lagrangians. The simplest example is furnished by 
the Jacobian determinant of an orientation preserving homeomorphism h G 
^i'"(X,Y) 

(1.50) [ J{x,h)dx = |Y|. 



One might ask which expressions of the type E{x, h, Dh) enjoy identities 
such as this? We call them free Lagrangians. Such a notion lies fairly 
deep in the topology of the mappings /i : X — t- Y. For example, using the 
topological degree, we find that the differential expression 

(1.51) E{x,h,Dh)dx = ^ ^^^^ 

is a free Lagrangian within the class of orientation preserving homeomor- 
phisms /i : A — )■ A* between anmili. Indeed, we have the desired identity 

(1.52) E{x, h, Dh) dx = Mod A for every /i G ^(X,Y) 
Ja 

A peculiarity of this example is further emphasized by the fact that the 
target annulus A* does not even enter into this identity. 

Like in the theory of polyconvex energies, the minimization of an energy 
integral whose integrand is a convex function of a number of free Lagrangians 
poses no challenge; Jensen's inequality usually gives the desired sharp lower 
bounds. Unfortunately, our integrand || Dh \\ " cannot be expressed as a 
convex function of free Lagrangians; though it is a convex function of the 
usual null Lagrangians. 



1.15. Uniqueness 

To reach the uniqueness conclusions, we first show that any extremal 
mapping /i : A — >■ A* satisfies the following system of nonlinear PDEs 

(1.53) D*h{x) Dh{x) = G{x, h) 

where the Cauchy-Green tensor G actually depends only on two variables 
I a; I and \h\. As this system is overdetermined it comes as no surprise that 
any two solutions h°{x) and h{x) for which |/i°(a)| = |/;.(a)|, at some a G A, 
must be equal modulo an orthogonal transformation T : — )• R", namely 
h{x) = Th°{x). Our proof exploits the classical computation of curvature 
of G in terms of its Christoff'el symbols. As regards the existence of such 
point a G A, we shall again rely on estimates of free Lagrangians. 



1.16. THE .ifi-THEORY OF INNER DISTORTION 
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1.16. The ^i-theory of inner distortion 

The conformal energy can naturahy be turned around so as to yield, the 
jSf^-integrabiUty of a distortion function of the inverse mapping [4]. 

We recall that a homeomorphism / : Y — t- X of Sobolev class W]^'^{Y, X) 
is said to have finite distortion if there exists a measurable function 1 ^ 
K{y) < oo such that 

(1-54) \Df{y)r^K{y)J{y,f) 

for almost every y G Y. Here \Df \ stands for the operator norm of Df. 
Using the Hilbert-Schmidt norm of matrices the outer distortion function of 
/ takes the form 

(1-55) K^(y,/) = i?Mi: 

^2 J(y,/) 

if J(y, /) > and we set ^^^(y, /) = 1, otherwise. There are many more dis- 
tortion functions of great importance in Geometric Function Theory. Among 
them are the inner distortion functions. 

4 



(1-56) K^iyJ) 
and 

(1-57) KJy, /) = 

^ ^ '^"-^^ ntdetL>«/(y) 

at the points where J{y, /) > 0. Otherwise, we set 

(1-58) K,iy,f) = K^{y,f) = l 

In recent years there has been substantial interest in mappings with inte- 
grable distortion [2, 4, 20, 23, 32]. Suppose now that / G l^^^f "^(Y, M"). 
The following identities connect the n-harmonic integrals with the theory of 
mappings with integrable inner distortion 

(1.59) / K^{y,f)dy= [ \Dh{x)\^ dx 

JY JX 

and 

(1.60) nt / K,(y,/)dy= / \\Dhix)\rdx 

where h denotes the inverse of /, [12, 19, 20]. These identities imply that 
the inverse of a mapping of integrable distortion always lies in the Sobolev 
space h G ^^'"(X, Y). Similarly, the ^^(A*, d;u)-integral means of K^iy, f) 
and (y, /) with respect to the dimensionless weight d/i = \y\~"' dy are the 
energy functionals 

^ ' Ja- W " Ja \Hx)\ 



■ dx 



^Here we use the operator norm of the cofactor matrix -D'/. 
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and 

(1.62) n2 / dy = / " ,, ,\; " dx 

Returning to (1.60), we infer from minimal n-harmonic mappings that 
Theorem 1.16. Under the Nitsche bounds 

(1.63) Kt(ModA) ^ ModA* ^ K^(ModA) 

the Sf^{Y)-norm of the inner distortion Kj (y, /) assumes its minimum value 
on a mapping f : A* A whose is inverse is a radial n-harmonic mapping 
h° : A A* . Such an extremal mapping f is unique up to a conformal 
automorphism of A. 

Theorem 1.17. // the domain annulus A* is too thin relative to the 
target annulus A; precisely, under the condition 

(1.64) ModA* < K|(ModA) -below the Nitsche bound. 

then the infimum of the ^^{A*)-norm of Kj{y, f) is not attained among 
homeomorphisms f : A* A. 

Nevertheless, we were able to find the infimum of the ^^-norms and 
identify the minimizing sequences. The weak BV-limits of such sequences 
and the underlying concept of their distortion (to be defined) are worth 
carrying out. 

Conclusion 

At the first glance the problems we study here may appear to be merely 
technical. However, their solutions require truly innovative approaches with 
surprising outcomes. For instance, the case n ^ 4 is different than one might 
a priori expect; upper bounds for the modulus of A* are necessary in order 
to ensure radial symmetry of the minimizers. The underlying technique of 
integration of various nonlinear differential forms is interesting in its own 
right. The free Lagrangians play a pivotal role, like null Lagrangians did 
play in the polyconvex calculus of variations. The entire subject grew out of 
fundamental questions of Quasiconformal Geometry about mappings of in- 
tegrable distortion. The paper also embraces a number of important aspects 
of the Calculus of Variation. 

We present this little theory here in two parts. 



Part 1 



Principal Radial n-Harmonics 



CHAPTER 2 



Nonexistence of n-harmonic homeomorphisms 



That nice smooth domains X and Y, such as annuh, may not admit a 
homeomorphism h : X Y of smahest conformal energy is a sequel of even 
more general observation. Let .4 be a topological annulus in R", we consider 
all possible n-harmonic homeomorphisms h: A ^ M"^. Unconcerned about 
the energy (finite or infinite) we address the question; which topological 
annuli A* can be obtained as images of A under such mappings? In general, 
this problem lies beyond our methods. Even in dimension n = 2 and when 
A* is a circular annulus (the Nitschc conjecture) the answer to this question 
required rather sophisticated ideas [26]. Nevertheless, one may roughly say 
that A* cannot be very thin. A specific instance is as follows. Let S be a 
convex {n — l)-hypersurface in R", given by 

S = {yGR":F(y) = 0}, 

where F is a "^^-smooth function in a neighborhood of S, such that VF ^ 
on S. We assume that the Hessian matrix 

" d^F 1 

— — -— is positive definite on S. 

dyidyij] 

Consider the e- vicinity of S; that is, = {y - \F{y)\ < e}. 

Proposition 2.1. If e is small enough then there is no n-harmonic 
homeomorphism h: A '-^ Vg such that S C h{A). 

Proof. We shall not give any explicit bound for e, but instead we argue 
by contradiction. Suppose that for every positive integer £, thus e = l/£, 
there exists an n-harmonic homeomorphism 

he = {hl,...,h'l): A°-^he{A), such that S c /i^(^) C V„ e = j. 

We appeal to the ^^'"-regularity theory of n-harmonic mappings [13, 36, 
47, 48, 49]. Accordingly, hi are uniformly bounded in "^^'"-norm on 
every compact subset of A. Moreover, there is a subsequence, again de- 
noted by {hi}, which converges together with the first order derivatives to 
a mapping h: A ^ S uniformly on compact subsets. The limit mapping 
h = (/i^, . . . , /i") still satisfies the n-harmonic equation. We aim to show 
that Dh = in ^, meaning that h is constant. Recall that Dh: A — > R"^" 
is (Holder) continuous. The computation below is certainly valid in the open 



21 



22 
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region where Dh ^ 0, because h is ^""-smooth in such region. We have 

F{h{x)) = for aU x e A. 
Applying partial differentiation and chain rule yields 



(2.1) 



E 



dF dh' 



for iy = 1,2, 



,n. 



^ % dxi, 

On the other hand, the n-harmonic equation takes the form 

n-2 dh' 



j/=i 

or, equivalently 



dxi, 



for z = 1, 2, . . . , n, where A = || Dh{x) 



(2-3) E 



(n - 2)A^ 



n— 3 



a/i^ . .„_2 9^^' 



dy, 



+ A 



dxiidx,. 



i 



,n. 



We multiply these equations by sum them up with respect to i, and use 
the identity (2.1) to obtain 



(2.4) 



,n-2 



^yi dxj^dxi, 

Next we diflFerentiate (2.1) with respect to u and sum up the equations, 



(2.5) 



EE 



dF d'^h' 



+E 



d'^F dh^ dh' 



i^it=i \ dy^ dxydxu ^ dyidyj dxy dx^ 



Here the double sum for the first term in (2.5) vanishes due to (2.4), so we 
have 



EE 



d'^F dh' dW 



dyidyj dx„ dx„ 



= 0. 



Since the Hessian matrix of F is positive definite, this equation yields ^ = 
for alH, v = 1,2, . . . , n, as desired. 

To reach a contradiction we look more closely at the homeomorphisms 
h^: A h£(A). Choose and fix an (n — l)-dimensional hypersurface T, C A 
(topological sphere) which separates the boundary components of A. Its 
images /i^(S) C he{A) separate the boundary components of he{A). Thus, 
in particular, 

< liminf [diam/i^(E)] = diam/i(E) =0 



a clear contradiction. 



□ 
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It is not so clear, however, how thick A* should be to ensure existence 
of n-harmonic homeomorphisms h: A A*. In dimension n = 2 the 
condition Mod^* ^ Mod^ is sufficient [24] and is not far from being 
necessary [27]. Questions of existence of harmonic diffeomorphisms between 
surfaces are treated in [33] . 



CHAPTER 3 



Generalized n-harmonic mappings 

In the classical Dirichlct problem one asks for the energy minimal map- 
ping h:X-^W of the Sobolev class h e K + #;^'"(X,M'') whose bound- 
ary values are explicitly prescribed by means of a given mapping ho G 
#'i'"(X,M"). Wc refer to the works of C. B. Morrcy [38, 39], where a 
systematic study of variational methods in the theory of harmonic integrals 
in vectorial case originated, see also [40]. The variation h h + erj, in 
which r) G R") and e ^ 0, leads to the integral form of the familiar 

n-harmonic system of equations 

(3.1) / {\\Dh II "-^D/i, Dt]) = 0, for every ry G '^^{X, M"). 
Jx 

Equivalently 

(3.2) Anh = Div ( [[ Dh \\ ""'^Dh) = 0, in the sense of distributions 

or, entry-wise, for h = {h^, . . . , /i") 

n 

(3.3) ^(||m[[--2/i^J^^=0, a = l,2,...,n. 

i=l 

In nonlinear elastostatics the matrix field 

S = Sh= II Dh{x) II "-2d/i(x) 

is known as Piola-Kirchoff tensor for the energy density function W{x) = 
II Dh{x) II This tensor represents the stress induced hy h:X^ M". Map- 
pings of Sobolev class y^^^"(X,M") that satisfy the equation Div<S^ = 
are called the equilibrium solutions. Equilibrium solutions in a given class 
ho + #i^'"'(X, M") represent unique minimizers within this class. The situa- 
tion is dramatically different if we allow h to slip freely along the boundaries. 
The inner variation works well in this case. This is simply a change of the 
independent variable; h^ = ho r]^, where r/^: X X are "^"^-smooth auto- 
morphisms of X onto itself, depending smoothly on a parameter e ~ where 
rjo = id: X X. Let us take on, as a initial step, the inner variation of the 
form 

(3.4) rieix) = x + eri{x), rj G '^^(X,W). 
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We compute 

Dhe{x) = Dh{x + er)){I + eDri) 
II Dh,{x) r = II Z^/i r + ne ( || Dh \\ ''-'^D*h ■ Dh , Drj) + o(e). 

Be aware that in this equation Dh is evaluated at the point y = x + e ri{x) G 
X. Integration with respect to x-variable yields a formula for the energy of 
he, 

Sf^^ = I [ II L>/i II " + ne{ II Dh \\ ''-^D*h ■ Dh , Drf)] dx + o(e). 

We now make the substitution y = x + er]{x) for which the following trans- 
portation rules apply: x = y — €r]{y) + o(e), r]{x) = r]{y) + o(l) and the 
change of volume element da; = [1 — eTr Dr]{y)] dy + o(e). The equilibrium 
equation for the inner variation is obtained from ^^he =0 at e = 0, 

(3.5) / {\\Dh II ■Dh--\\Dh II , Dt]) dy = 
or, by means of distributions 

(3.6) Div ^ II Dh II "-2r>*/i -Dh-^jDh II "7^ = 0. 
Now we introduce the divergence free tensor 

A= \\Dh\\''-'^D*h- Dh--\\Dh\\''I = ( C --TiC\Tr'^C 



n \ n 



where we recall the right Cauchy-Green tensor C = C{x) = D*h ■ Dh. The 
name generalized n-harmonic equation will be given to (3.6) because of the 
following: 

Lemma 3.1. Every n-harmonic mapping h G >r;^;"(X,M") solves the 
generalized n-harmonic equation (3.6). 

Proof. We consider again the perturbed mappings /ie(x) = h{x+e r]{x)) 
which coincide with h outside a subdomain U compactly contained in X. 
Applying the integral form (3.1) of the n-harmonic equation, but with the 
test function h — h^ & Wo '"(1LJ,M") in place of r], ^ we estimate the energy 
of h over U as follows 

<4= / \\Dh\\'^= I \\Dhr-^{Dh, Dh) 
Jv Jv 

n-l 1 

= J \\Dhr-^{Dh, Dh,) ^(^j \\Dhr^ " \\Dh,ry . 

Hence 

Sh ^ S'he 1 with equality at e = 0. 
This means = at e = 0. Equivalently, h satisfies the generalized 

n-harmonic equation. □ 



^This is justified because '^o°°(U,R") is dense in >ro^'"(U,R"). 
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In dimension n = 2, the generalized harmonic equation reduces to 

(3.7) Div (^D*hDh -^\\Dh II = 0. 

This equation is known as Hopf-Laplace equation [11]. In complex notation 
it takes the form 

d 

(3.8) — {hji^ =0, z = xi+ 1X2. 
If, by some reason, h G "^^(X, C) then (3.8) reads as 

(3.9) Jiz, h)Ah = 0, where A = 4——. 

ozoz 

There exist diverse non-harmonic solutions to (3.9), many of them of great 
interest in the theory of minimal surfaces and some with potential applica- 
tions to nonlinear elasticity (elastic plates), see [11]. 

Equally in higher dimensions, the n-harmonic mappings of Sobolev class 
#^^^"(X, M") are only particular solutions to the generalized n-Laplacian. 

3.1. Solutions to the generalized n-harmonic equation that are 

not n-harmonic 



h{x) = H{jx\)-r—r, where H = H{t) is absolutely continuous 



Let us take on stage a radial mapping 

X 

\x\ 

see Chapter 4 for extensive treatment of such mappings. We find that 
(3.10) 

A = \\Dh \\ (D*h-Dh-- II Dh II 2/ 

n 



(n-1)^ //2 + ^ 



V / l^l" \ \x\^ n J 



It is shown in Chapter 5, that if ^ is a "^^-smooth n-harmonic mapping then 
H = H(t) must satisfy the characteristic equation 



n-2 
2 



(3.11) (ij2^^!^J " . _ |^|2^2j ^ 



const. 



This also confirms that h satisfies the generalized n-harmonic equation be- 
cause of the identity 



DivJ-f:^-i/)=0. 

x r n 



However, the hammering mapping h[x) = j^, corresponding to i^(t) = 1, 
solves equation (3.6) but is not n-harmonic. If we paste this hammering 
mapping with a compatible smooth radial n-harmonic solution (the critical 
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Nitsche mapping) there will emerge a ^^'^-smooth minimal deformation be- 
tween annuli. There are many more non n-harmonic solutions to (3.6), some 
seem to be unsatisfactory; for example those with Jacobian changing sign are 
forbidden in elasticity theory because of the principle of non-interpenetration 
of matter. 

3.2. Slipping along the boundaries 

Let us return to the inner variation in (3.4), he{x) = h{x + erj^x)), but 
this time with ij not necessarily having compact support. We assume that X 
is a "^^-smooth domain and rj G "^^(X). The fact that h is allowed to freely 
slip along the boundaries amounts to saying that the vector field 77 = r]{x) 
is tangent to 9X at every point x G 9X. The integral form of the resulting 
variational equation is the same as in (3.5). Integration by parts (Green's 
formula) will produce no integral over X, because of (3.6); there will remain 
only boundary integrals. Precisely, a general formula we are referring to is: 

0= / {A,Dr])dx = - / (DivA, 77) dx / (A(x)n(x), ?7(x)) dcr(x) 

Jx Jx JdX 

where n{x) is the outer unit vector field to 9X and da{x) stands for the 
surface measure. This is justified under an appropriate assumption on the 
degree of integrability of A, DivA, rj and D-q. Nevertheless in our case, we 
obtain 

(3.12) / (( II Dh II ■Dh--\\Dh II "/)n(x) , rj{x)) da{x) = 
JdX n 

Since rj{x) can be any tangent field the equation (3.12) is possible if and 
only if the vector field 

^ II Dh II ''-^D*h -Dh-^WDh II '^/^ n(x) 

is orthogonal to dX. In other words, at each x G dX the linear mapping 
II Dh II ^~'^D*h- Dh — i II Dh \ "/ takes the (one-dimensional) space of normal 
vectors into itself. Of course, the same holds for the mapping D*h-Dh. Since 
D*h ■ Dh is symmetric one can say, equivalently, that D*h ■ Dh preserves 
the tangent space. We just proved the following. 

Proposition 3.2. The equilibrium solution for mappings that are slip- 
ping along the boundaries satisfies, in addition to (3.6), the following con- 
dition 

(3.13) D*h- Dh: T^dX^T^dX 
equivalently, 

(3.14) D*h ■ Dh : N^dX N^dX. 

where TxdX and TSixdX designate the tangent and normal spaces at x E dX. 
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In dimension n = 2 this amounts to saying that the Hopf quadratic 
differential hzhz dz^ is real along the boundary components of X. In an 
annulus K = {z: r < \z\ < R} every Hopf differential which is real on dK 
takes the form 

— c 

hzhz = c is a real constant 

z^ 

see [11]. One of the solutions is the mapping we are already encountered 
in (1.20) 

{jf| ;^ < l-^l ^ 1 ) hammering part 

i + =) 1 ^ l^l ^ , harmonic part 

In higher dimensions the radial mapping h{x) = H[\x\^-^ in the annulus 
A = {x: r < \x\ < R} also complies with the boundary condition (3.13). 
Indeed, the Cauchy-Green tensor takes the form 

Fr V \^\ J fI 
The normal vector field at is n = Xi{x) = x, and we see that 

[D*h{x) ■ Dh{x)] n = An, A = [i^d^DJ^. 

For any tangent vector t = t{x), we have {x (g) x)t = so 

[D*h{x) ■ Dh{x)] t = ^t. 

The corresponding singular values of the right Cauchy-Green tensor are 
H{x) and — principal stretches. 

Before moving to a study of n-harmonic mappings between spherical 
annuli, let us fulfill the promise of proving that weak /^^'"-limits of home- 
omorphisms are indeed monotone. 

3.3. Proof of Theorem 1.7 

We are dealing with a sequence of homeomorphisms hj : X Y, j = 
1, 2, . . . converging c-uniformly to a mapping /i: X — >■ Y such that 

dist(/ij(x), (?Y) ^ ri{x) II Dhj \\ ^n(x)- 

see (1.16). Passing to the limit we also have 

dist(^(x), 5Y) ^ r){x) liminf || Dhj \\ ifn(x)- 

Therefore, for every e > there is 5 > such that 

(3.15) dist{hk{x),dY) < e whenever dist(x,5X) < S 

for every A; = 1, 2, . . . and 

dist(h(x), 5Y) < e whenever dist(a;, 5X) < S. 
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This shows that h^^(K) C X is compact whenever K C Y is compact. To 
show that h is also monotone we argue as follows. Consider a continuum 

K C Y and assume, to the contrary, that /i~^(]fC) = A U B, where A and B 
are nonempty disjoint compact subsets of X. For every integer i such that 

l<eiiidist(K,5Y) 

one may consider a neighborhood of IK, 

Ke = {ye Y: dist(y,K) < l/£}. 

This is an open connected set in Y. The preimage h~^{K£) C X is a neigh- 
borhood of the compact set A U B C X. Since hj h c-uniformly then for 
sufficiently large j, say j = je, we have 

The sets h'^iKi) are connected. Therefore to every £ ^ j there corresponds 
a point X£ G C X such that 

(3.16) dist(a;£,A) = dist(x^,B). 

We also note that 

dist{hj^{xe),dY) ^ dist(]K^,aY) ^ 2e - j > e. 

Hence, in view of condition (3.15), we have 

dist{xi,dX) ^ 5, £ = 1,2,... 

We then choose a subsequence, again denoted by x^, converging to some 
point a; G X. Passing to the limit in (3.16) yields 

dist(x, A) = dist(x,B), thusx^AuB. 

On the other hand, since hj converge to h uniformly on compact subsets, 
we may pass to the limit with the sequence hj^{xi) — ^ h{x). Finally, since 
hj^{xi) G IK^ and PlK^ = K, we conclude that h{x) G K, meaning that 
X G h~^{K) = A U B. This contradiction proves Theorem 1.7. 



CHAPTER 4 



Notation 



This chapter is designed to describe basic geometric objects and concepts 
that will be used throughout this paper. The primary domains here are 
annuli in M", also called spherical rings. They are subsets of the punctured 
Euclidean space, 

=R"\{0} = !^X = {xi,...,Xn); \x\ = yjxl + ---+xly^0^ , 71^2 

4.1. Annuli and their modulus 

There are four types of annuli, generally denoted by A C Mq 
(4.1) 

A(r, i?) = {x G R" ; r <\x\< R] , where ^ r < ^ oo 
^ _ ^ A[r, R) = {x ; r ^\x\<R} , where < r < ^ oo 
~ I A(r, i?] = G R"- ; r < |x| ^ i?} , where ^ r < i? < oo 
A[r, i?] = {x G R" ; r ^ |a:| ^ i?} , where < r < 7? < oo 

The boundary of each annulus consists of two components (the inner and 
the outer sphere) except for the degenerate case A[r, r] which reduces to the 
{n — l)-dimcnsional sphere denoted by §,"^^. 

The conformal modulus of an annulus is defined by 

R f dx 

(4.2) ModA = a;„_ilog- = -— g[0,oo] 

r Ja 

Hereafter, Un-i denotes the surface area of the unit sphere S*^"^ C R". 

(4.3) COn-l = 



''^'""'^ for n = 2A; + 1 



l-3---(2fc-l) ' 

4.2. Pol£ir coordinates in R" 

Polar coordinates are natural to use when working with annuli. We 
associate with any point x G R" a pair of polar coordinates 

(4.4) {r,oj) G R+ X S"-i ~ RJ^ 

where r = \x\ is referred to as the radial distance and a; = A as the 

I I \x\ 

spherical coordinate of x . Obviously x = rui and the volume element in 
polar coordinates reads as da; = r'^^^drdo; 
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4.3. Spherical coordinates, latitude and longitude 

A position of a point w G S""^ is determined by its latitude and longi- 
tude. These are convenient coordinates for describing deformations of the 
unit sphere. We refer to the points a;+ = (0, ...,0, 1) and = (0, ...,0,— 1) 
as north and south poles, respectively. The unit sphere with poles removed 
will be denoted by S!^~^. The equatorial sphere S""^ C S"~^ is given by the 
equations + ... + a;^_i = 1 and Xn = 0. To every point u G S^~^ there 
correspond meridian coordinates 

(4.5) {e,s) G (0,7r) X S"-2 ~ S^-^ 

Here s lies in the equatorial sphere while 6 is the distance south of the north 

pole measured in degrees. The rectangular coordinates of a; = (xi, X2, Xn) 
are uniquely recovered from s and 9 by the rule 

(4.6) a; = (cos6',s-sin6'), < 9 < tt 

The (n — l)-surface area of S"~^ is expressed in terms of the meridian angle 
as 

(4.7) do; = (sin^)"-2 d^ds , / du = ojn-i- 

where ds stands for the (n — 2)-surface area of the equatorial sphere. There- 
fore, 

(4.8) rsin-2^de = ^ 

Jo 

We turn now to the basic examples of mappings between annuli. 



4.4. Radial stretching 

These are the mappings which change the radial distance but leave the 
spherical coordinate fixed. Let s^C [r,R], ^ r ^ R < oo, denote the space 
of absolutely continuous functions on a closed interval [r, R] . Associated 
to every H G AC [r, R] is a mapping /t : A — > A* of the annulus A = 
{x ; r ^ \x\ ^ R} onto an annulus A* = {y ; r* ^ \y\ ^ R*}, defined by 
the rule 

(4.9) h{x) = H{r) u , where r = \x\ and w = -j^ G S"""^ 

\x\ 

Here we have 

(4.10) r* = min \H{t)\ ^ max \H{t)\ = R^ 

We call h a radial stretching and refer to H : [r, R] — t- [r^, R^] as its normal 
strain. Note that a composition of two such mappings is again a radial 
stretching; absolute continuity is not lost. 



4.5. SPHERICAL MAPPINGS 
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4.5. Spherical mappings 

A generalization of radial stretching is the so-called spherical mapping 

(4.11) h{x) = H{\x\) $(^) , where $ : S""^ ^ S""^ 

\x\ 

Here, as in the case of radial stretchings, the normal strain function is as- 
sumed to be absolutely continuous. The spherical part <1> : S"^^ — ;> 
called tangential tension, is continuous and weakly diffcrcntiablc. In what 
follows we will be actually concerned with homeomorphisms $ : S"^^ — )■ 
gn-i ^]^g Sobolev class W^'"'{S'^'~^,S'^~^). Composition of two spherical 
mappings results in the composition of their normal strain functions and 
tangential tensions, respectively. 



CHAPTER 5 



Radial n- harmonics 



In this chapter we identify all solutions to the n-harmonic equation 

(5.1) div II Dh II ""-^Dh = of the form h(x) = H(\x\) ^ 

Such solutions, called radial n-harmonics, will be originally defined in the 
annulus 

(5.2) A = A(a,6) = {x; a<|x|<6}, 0<a<6<oo 

Concerning regularity, we assume that the strain function H belongs to the 
Sobolev space 

1 Tl — 1 

(5.3) <f(a,6)c^aa,fc), " = ^ 

We are going to give a clear account of how to solve the n-harmonic equation 
for the radial n-harmonics. But before, let us state in advance that the radial 
n-harmonics will actually extend as n-harmonics to the entire punctured 
space Rq. Even more, they will extend continuously to the Mobius space 
MP = M" U {oo}, one point compactification of M". 

Let Tin denote the class of all n-harmonics in M", the following elemen- 
tary transformations of variables x and h preserve this class. 

• Rescaling; 

h{x) £ Tin implies A h{kx) G Tin for every k > and A G M 

• Reflection; 

h{x) G T-Ln implies h (xjxj^^) G Hn 

The radial n-harmonics inA = {x; o<|a;|<6},0<a<6<oo are none 
other than the unique solutions of the Dirichlet problem 

' div II Dh II ""-"^Dh = for a < |x| < 6 

(5.4) < h{x) = ax for |a;| = a 

h{x) = fix for |a;| = h 

where a and ^ can be any real numbers, so H{a) = aa and H{h) = fib. It is 
understood here that h extends continuously to the closed annulus A [a, 6]. 
We shall distinguish four so-called principal n-harmonics in and use them 
to generate all radial n-harmonics via rescaling of the variables x and h. For 
this we need some computation. 
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5.1. The n-Laplacian for the strain function 

Suppose H G "^^(0,6), ^ a < 6 ^ 00. Thus the derivatives of the 
radial stretching 

(5.5) h{x) = H{\x\)^ 

\x\ 

exist up to order 2. The differential matrix of h can be computed as 

(5.6) DH^) = :5^Ml + . ^ 

\x\ \x\ \x\'^ 

Hereafter H stands for the derivative of H = H{t), a < t < b. Hence the 
Hilbert-Schmidt norm of the differential matrix is expressed by means of H 
as 

tt2 

(5.7) II Dh(x) II ^ = Tr (D*h ■ Dh) = ii"^ + (n - 1)^ 

|x|^ 

A lengthy computation leads to the following explicit formula for the n- 
Laplacian 

where t = \x\. Our subsequent analysis will be based on the elasticity 
function already mentioned at (1.30), 

(5.9) ^ = ^^ = ^- 

Let us make certain here that we shall never deal with the case in which 
both H and H vanish simultaneously, so will be well defined even when 
H vanishes. Now h is n-harmonic if and only if 

(5.10) (1 + rf') t^H = (1 - r?) [rf" + {n - 2)r] + n - l] ■ H 

where the case 7] = 00 should read as H = 0. In dimension n = 2 the 
n-harmonic equation reduces to the familiar linear Cauchy-Euler equation 
for H 

(5.11) t'^H + tH - H = 

It has two linearly independent solutions, denoted by 

(5.12) Ho{t)=t and Hoo{t) = ^ 

They generate all solutions by means of linear combinations. We pick out 
two particular solutions 

(5.13) H+it) = Ut + ]) and H^it) = ^ ' ^ 
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and refer to Ho, H^o, and as principal solutions. The reason for 
preferring these principal solutions to the ones in (5.12) is simply that we will 
be able to generate all radial n-harmonics, simply be rescaling the principal 
ones. Such a set of principal solutions will be particularly useful in higher 
dimensions where the n-Laplacian is nonlinear. For n = 2 the corresponding 
complex principal harmonics are 

(5.14) hoiz) = z, hooiz) = ^, 

(5.15) h+{z) = ^(^z + ^^ and h.{z) = ^ (^z - ^ 

The reader may wish to observe that ho and hoo map univalently onto 
itself, whereas h+,h- have a branch set (folding) along the unit 

circle. Precisely, covers the exterior of the unit disk twice while h- is 
a double cover of the entire space R^. In spite of nonlinearity the analo- 
gous ordinary differential equations in higher dimensions arc still effectively 
solvable, with some computational efforts. The key is the following identity 
which follows from (5.8). 

Lemma 5.1. Leth{x) = H{\x\) be a radial stretching of class 'Tf'^{A,W'), 
where A = {x; a < |x| < b}. Then, with the notation \x\ = t , we have 

""^■div \\Dhr-^Dh = 



n-1 

rl TO — 2 

(5.16) + ^ {H^-t'H^)} 

5.2. The principal solutions 

The main idea is to eliminate one derivative in the second order n- 
harmonic equation. This will lead us to a problem of solving a first order 
ODE with one free parameter. For nonlinear equations such a procedure is 
rather tricky. Fortunately, there is a very satisfactory realization of this idea 
due to the identity (5.16). We wish to express the n-harmonic equation for 
h{x) = H{\x\) 1^ in the form 

(5.17) CH = F{t, H, H) = const. 

By virtue of Lemma 5.1, the nonlinear differential operator C takes the 
following explicit form 



(5.18) CH 



H^ + 

n — 1 



n-2 

2 



'h^ - t^H^^ = const. 



^In general, guessing the identity such as (5.18) takes some efforts; though verifying 
it poses no challenge. 
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Certainly, all radial n-harmonics can be found by solving the so-called char- 
acteristic equation 

(5.19) CH = const. 

However, for the converse statement, caution must be exercised because 
this equation only implies that H div \\ Dh\\'^~'^Dh = 0. For example, the 
constant function H{t) = 1 satisfies CH = 1, but the corresponding radial 
mapping h{x) = |^ is not n-harmonic. ^ More generally, there are many 
more solutions to the equation CH = 1 which are constant on a subinterval. 
For example, one can paste with a suitable n-harmonic map on the rest 

of the interval. Such solutions, however, cannot be "^^-smooth; '^^'"'^-regular 
at best. They, nevertheless solve the generalized n-harmonic equation (3.6). 




Figure 1 . A solution which is constant on an interval. 

Solutions H = H{t) to the characteristic equation CH = c, that are 
constant in a proper subinterval like in the above figure, cannot be even 
■^i.i.j-egular if c ^ 0. 

Let us summarize our findings as: 

Definition 5.2. The term principal solution pertains to each of the 
following four functions of class ^^(0, oo) which solve the equation CH = 

constant: 



(5.20) 

(5.21) 

(5.22) 
(5.23) 



Ho{t)=t, 

Hoc{t) = J 

CH+ = 1, 
CH^ = -1, 



CH^ = 

CH^ = 

H+{1) = 1 
iJ_(l) = 



We shall solve these Cauchy problems and show that both H+ and 
are actually "^""-smooth. 



This noninjective solution will play an important role in the sequel. 



5.3. THE ELASTICITY FUNCTION 
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Figure 2. The principal solutions for n = 2. Similar ones 
for n ^ 3 are illustrated in Figures 7 and 10. 

The name "principal solutions" is given to ifo, H^, and be- 
cause they generate all radial n-harmonics. We denote the corresponding 
radial mapping by ho, h^o, ^+ and h-, and call them principal n-harmonics. 
Precisely, we have 

Proposition 5.3. Every radial n-harmonic mapping in the annulus A = 

{x; a < |x| < 5} takes the form 

(5.24) g{x) = A • h{kx) , XgR, k>0 
where h G ^°°{Wl), is one of the four principal n-harmonics. 

It should be observed, as a corollary, that radial ra-harmonics are 
smooth in the entire space . 

5.3. The elasticity function 

We shall distinguish four classes of the radial ra-harmonics h{x) = 
The concept of so-called conformal elasticity underlines this distinction. In 
what follows it will never be the case that both H and H vanish simulta- 
neously. Consider any function H G ^^(r, i?), < r < R < oo, such that 
H"^ H"^ > 0. The elasticity function 

(5.25) r,{t) = r,,{t) = ^-§^ 

is continuous on (r, R) . In general, it assumes values in the extended real 
line R = M U {-oo, oo}. We recall that h = H{\x\) ^ is: 
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• Conformally contracting, if \ri{t)\ < 1 

• Conformally expanding, if \r){t)\ > 1 

• Conformal, if \r){t)\ = 1 

Observe that the equation r]{t) = 1 yields H{t) = X ■ t, while r]{t) = —1 
gives H{t) = J. These are special cases of the power functions H(t) = At", 
having constant elasticity r/(t) = a. The above three classes are invariant 
under rescaling and inversion. Precisely, we have the formulas: 



(5.26) 



\r]jj{kt) = r]p{t) where F{t) = H{kt) 

I {l/t) = -ri^ (t) where F{t) = H{l/t) 



The elasticity function tells us something about the infinitesimal relative rate 
of change of the modulus of an annulus under the deformation h = i? (|x|) 
Precisely, we have foi r < t < R 

_ ModA*{t + e)^ModA*{t) 
(5.27) [t) - lim ^ _ ^^^^ 

where A{t) = A{r,t) and A*{t) = h[A{t)]. 

We shall be concerned with ^^-solutions to the characteristic equation 



(5.28) £H 



n — 1 



This equation is invariant under rescaling and inversion of the independent 
variable t. Accordingly, all ^^-solutions fall into three categories: 

• H is contracting if c > 0, equivalently \rj{t)\ < 1 

• H is expanding if c < 0, equivalently |??(t)| > 1 

• H is conformal if c = 0, equivalently \ri(t)\ = 1, 
Furthermore, in the conformal case 

( h is preserving the order of boundary components if r]{t) = 1 , H(t) = Xt 
\ h is reversing the order of boundary components if ri{t) = —1, H{t) = j 



This is easily seen by writing (5.28) as 

n — 2 

We then distinguish four classes of the radial n-harmonics: 



(5.30) 
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where 



= 


1"= 


H{\x\)-^ 


U- = 




H{\x\)- 


Ho — 


[h; 


h{x) = 


^oo = 


ih- 


h{x) = 



X 



x\ 



A • x 



A G 



We strongly emphasize again that not every solution to the equation jOH = 
c defines a radial n-harmonic map, due to lack of regularity. There are 
Lipschitz solutions to CH = lacking ^^-regularity 




Figure 3. Lipschitz solution to £H = 0. 



Such solutions are not particularly desirable. Even ^ -solutions to the 
equations CH = 1 need not give radial n-harmonics. For example, if we 
paste Hj^{t) = 1(^ + 7) l^i<oo with function identically equal to 1 for 
< t ^ 1 then the resulting function will become a ^^'^-solution to the 
equation CH = 1 on (0,oo). 




1 



Figure 4. A ^^-solution lo CH = 1. 
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The radial stretching obtained in this way fails to be n-harmonic on 
the punctured disk A (0,1). Thus we see various degenerations of 
sohitions, suggesting that we should restrict ourselves to the solutions of 

class <^2(o,oo).3 

The two principal solutions iJ+ and i7_ are particularly interesting and 
important to the forthcoming results, so we devote next two sections for a 
detailed treatment. 



5.4. The principal solution (conformal contraction) 

We find i7+ as a "^^-solution to the Cauchy problem 



(5.31) 



H 



n 



(h^ - fH'^^ = 1 , H{1) = 1 



The general idea behind our method of solving this equation is as follows. 
Differentiating (5.31) yields a first order equation for the elasticity function. 



(5.32) 



(1 + 77^)?) 



(1 - r/2) {rj'^ + n-1) t ' 



m\ < 1 



This can also be directly derived from (5.8). With the equation (5.32) at 
hand we now proceed to the explicit computation. Consider the function 
r+ = r_|_(s) defined for — 1 < s < 1 by the rule 



(5.33) 



(1 + r2) dT 



')(t2 + n- 1) 



T^f exp 



tan"^ ' 



_n^/n—l vn— 1_ 

Obviously r+ is strictly increasing from zero to infinity, thus invertible. 



f r. 




Figure 5. r+ is strictly increasing and "^""-smooth. 
'In the expanding case 'if^ solutions are always smooth. 
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Note the identities 



(5.34) 



r+(o) = 1 

r+(o) = 
r+(-s)r+(s) = i 

We then examine the inverse function defined for < t < oo, 

(5.35) u{t) = u+{t) = T-^{t), u(l) = 0, 'ii(l) = n-l 

The identity r+(— s) r_|-(s) = 1 translates as the antisymmetry rule for u 



(5.36) 



-u{t) , 



t > 



Implicit differentiation of the equation r+(u(t)) = t yields 

r+(u) 



(5.37) 



■uf+(tt) = 1 



tu{t) 



f+(u) 



We shall now introduce = H^{t), < t < oo, and prove that iJ+ is the 
principal solution, namely 



(5.38) 



1 + 



n — 1 



n 2 



(l - U^) " ^ 1 



Prom this formula it follows that u is the elasticity function of Indeed, 



(5.39) ri = ri+it) = 



H 



+ 



= u - tu - — — = u{t) 



(1 - u^){n - 1 + «2) 



'r+(n) 




r] = r]^{t) 



Figure 6. The elasticity function of 
Now, the equation CH^ = 1 can easily be verified 



CH. 



+ 



(5.40) 



\H. 



+ \ 



1 + 



1 + 



n 



,2 1 



n-1 



n-1 



(l-u^).l 
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This last step is immediate from (5.38). It is worth noting that G 
^°°(0, oo) and -ff+(l) = 0. Antisymmetry rule at (5.36) results in a sym- 
metry rule for 



(5.41) 




= H+ {t) 



Furthermore, it follows from (5.10) that H+{t) > 0. Thus is strictly 
convex, see Figure 7. 




Figure 7. The principal solution H+. 



Let us take a look at the behavior of i?+ near zero and infinity. It has 
an asymptote at infinity with the slope 



n-2 
1 \ 2n 1 

(5.42) e = e+ = ( 1 - - j 4-n exp 



n-2 



n-\/n — 1 



tan 



1 



Indeed, under the notation above, we have the following equation for t 



1 + 



n- 1 



Note that A = A{s) is "^""-smooth on (— l,oo). It vanishes at s = 1. 
Application of L'Hospital rule shows that 



lim [i7"(t)-e"t"]=y4'(l) >0, 



(5.43) 

whence the inequality 
(5.44) 

In particular by (5.41) we conclude that 
(5.45) 



(}<H{t)-et ^ 



rv 



Qn-l 



lim = e = MmtHit) 

t-^oo t 
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5.5. The principal solution i7_ (conformal expansion) 

The principal solution H = iJ_ G ^°°(0, oo) is obtained by solving the 
following Cauchy problem for the characteristic equation 



(5.46) 



n ■ 



1 



n-2 
2 



[h'^ - fH'^^ = -1 , H{1) = 







Obviously H cannot vanish. For an explicit computation we introduce the 
following function 



(5.47) r_(s) = 



l + s 



■ exp 



n-2 



n 



tan ^ (^s^/rl^^T) 



where — 1 < s < 1. Its logarithmic derivative is computed as 



(5.48) 



r_(5) 



> 



r_(s) (1 - s2) [1 + (n - l)s2] 
This shows that r_ : (—1,1) — >• (0, oo) is strictly increasing. Consequently, 



we consider its inverse function 



(5.49) 



u{t) = u-{t) = r-^(t). 



< t < oo 



The graphs of r_ and u_ are very similar to those of r_|- and u^, though 
they are not the same. However, this time the same identity r(— s) r(s) = 1 
yields different transformation rule for u, 



(5.50) 



u I - 



Mi) , 



< i < oo 



It is at this point that the analogy with ends. We define the principal 
solution iJ_ as 



(5.51) 



H{t) = H_{t) = u (u' + 



1 



n 



1 



n 2 



Note that i7_ changes sign at t = 1, exactly where u vanishes. Basic proper- 
ties of H-{t) can be derived directly from (5.51). Let us state some of them. 
We observe the first difference between iJ+ and iJ_ in the antisymmetry 
formula 

^ ' =-H-{t) 



(5.52) 



which is immediate from (5.50). Next for t ^ 1 we compute the elasticity 
function for 



rj = rj_ (t) 



(5.53) 



tH-{t) _ {l + u^)tu 

H^{t) ~ u(l-n2)[l + (n-l)n2] 
til f_(u) _ 1 

u r_('u) u 
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' V = V- (t) 

Figure 8. The graph of the elasticity function for H^. 



Now the characteristic equation for H = H-{t) is easily verified 



(5.54) 



1 + 



n — 1 



u 



1 



(1 - V') 
1 



n — 1 



71-2 

2 



Our analysis of convexity properties of is based on the equation (5.10), 
which we write as 

(5.55) (1 + u^)t^ H = {u-l) [{n - l)u^ + (n - 2)u + l] — 

for alH 7^ 1. Since the quotient ^ is positive we see that 

(5.56) sgnH- = — sgn [(n — l)u^ + (n — 2)u + l] 

The key observation is that for n = 2, 3, 4, 5, 6 the polynomial P = Pn{u) = 
(n — l)u'^ + (n — 2)n + 1 has no roots, meaning that iJ_ is concave. Rather 
unexpectedly, is no longer concave when n ^ 7, because P„ has two 
roots in the interval u G (—1,0). 



f P 




P = Pn{u) 




Pn{u) 



Figure 9. Pn is positive if 2 ^ n ^ 6 and has two roots if 
n>7. 
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As a consequence, the graphs of H- exhibit two inflection points when 
n ^ 7. 



n = 2,3,4,5,6 




H = H_{t) 




2 inflection 
points 



Figure 10. The principal solution H_. 



Analysis of the asymptotic behavior of iJ_ is much the same as that for 
H^. We only state the results. The slope of the asymptote at infinity equals 



n-2 

1 4 n exp 

n I 



(5.57) e = 0- 

and we have 
(5.58) 

In particular, for H = H_{t) it holds 
(5.59) 



2-n 



<et- H(t) ^ 



n\/n — 1 
C 



tan ^ ^Jn — \ 



lim = e = - XmitRit) 

t-^oo t t->0 



5.6. The boundary value problem for radial n-harmonics 

Given a, /? G M one may wish to look for functions H G "^^(a, 6)n^[a, 6], 
< a < 6 < oo, which satisfy the conditions 



(5.60) 



£H = c 

H{a) = a 



Here the constant c is also viewed as unknown, otherwise the system (5.60) 
would be overdetermined (ill-posed). Let us begin with the easy case (3 = 0. 
The solution is given by the formula 



(5.61) 



H(t) = \H_{kt), 



b 



A 



a 
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From now on we assume that /3 / 0. There are five cases to consider. 
Case 1. Suppose that 

(5.62) — < ^ < 5 

The solution wiU be found in the form H(t) = XH-{kt). One has to show 
that there exist real number A and a positive number k such that 



(5.63) 



\XH_{ka) =a 
\xH_{kb) = P 

We eliminate A by dividing the equations 

Note that hm Q{k) = — oo. On the other hand using (5.59) we see that 



I H-{ka) \ ^ 



(5.65) hm Q(k) = - hm „ ,,,, , _ _ 

^ ' fc^oo^ ^ b \k^oo H-jkb) / 5 5 

\ kb / 

By Mean Value Theorem there exists cc > ko > ^ such that 

H-{koa) a 



(5.66) Q{ko) 



Finally, conditions (5.63) are satisfied with A = -jf^^^- 

The next case is obvious; a linear function is a solution. 
Case 2. Suppose 

hi 

Then the solution is given by 

(5.68) H{t) = p = aH,{i^) 

Case 3. Suppose now that 

a a b 

We are looking for the solution in the form H(t) = XH^[kt). As before, one 
needs to find ko such that 

H+{koa) _ a 
H+{kob) - /? 
For this reason we consider the function 

With the aid of the asymptotic formula (5.45) we see that ^lim Q{k) = | 

and lim Qik) = -. By Mean Value Theorem the equation (5.70) holds for 
k^o " 
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some ko > 0. Then the parameter A can be chosen to yield the boundary 
conditions at (5.63). 

Next case is immediate. 
Case 4. Suppose 

The solution is given by the formula H{t) = ^ = a (|) . 
Case 5. Suppose finally that 

(.73) ^ > t 

We are looking for the solution in the form H(t) = XH-(kt), so we examine 
the function 

By virtue of formula (5.59) we see that lim Q{k) = +oo and lim Q{k) = -. 

This case is completed by again invoking Mean Value Theorem; there is 
ko < I such that 

H_{koa) _ a 

^^•^^^ Hlik^bj ~ P 

Then, we take A to ensure that XH-{kob) = p. 

In summing up this section, we note that the foregoing analysis, com- 
bined with the uniqueness of the Dirichlet problem for the n-harmonic equa- 
tion, reveals that in fact the principal n-harmonics generate all radial n- 
harmonics. As an unexpected bonus we deduce that every radial n-harmonic 
map in a ring domain extends n-harmonically to the entire punctured space 
M» = \ {0} and continuously to the Mobius space = R'* U {oo} . 



CHAPTER 6 



Vector calculus on annuli 



In this chapter we briefly review differential calculus on annuli. 



6.1. Radial and spherical derivatives 

In studying deformations of spherical rings one must consider the radial 
and spherical components of a differential map. Through every point x G M" 
there passes a sphere 



(6.1) Sr' = {yeR"; \y\=t}, t = \x\ 
Its tangent hyperplane at x is given by 

(6.2) T, = {e; {x,0=t^} 

Fix an orthonormal basis for T-^, say T = {T2,T^, ...,Tn}} Let O C M" be 
a domain and ^ : J7 ^ R" a mapping having first order partial derivatives 
defined at x G The radial (or normal) derivative of ^ at x is a vector 
defined by the rule 



h^{x) = Dh{x) ■ N 



^g_3') u _ 7^^,^■™^ tvt _ ^^^xi + ••• + 

where N = is called the radial (or normal) vector. In an exactly similar 
way we define the spherical (or tangential) derivatives 



(6.4) 



hTix) = [hT2,---,hTj , hr^ = Dh{x) ■ Tt, , u = 2,3,...,n 



This is an (n — l)-tuple of vectors in R", conveniently considered as column 
vectors of an n x (n — l)-matrix. Continuing in this fashion we view the dif- 
ferential of Dh as an n X n-matrix whose column vectors arc /?,n, hx.2 , /it„ ■ 
Of course, such a matrix representation of Dh depends on choice of the or- 
thonormal frame N, T2, T^j 



(6.5) 



Dh{x) 



/iN tlTo 



^It is not always possible to choose such a basis continuously depending on a; e R" 
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However, some differential expressions do not depend on the frame. For 
example, the Hilbert-Schmidt norm of the differential matrix, 

(6.6) II Dh\\^ = |/iN|' + \hT, f + ... + |/iTj' 



Also, the spherical component of || Dh \\ , which we define by the formula 

(6.7) iH = f:l"""'T ^(i'"-i---i'""0'^ 



is frame free. Equality holds if and only if I/1T2I = ... = |^T„|. Thus we have 

(6.8) \\Dhf = \h^f + {n-l)\hTf 

Another differential quantity of interest to us is the Jacobian determinant 

(6.9) J{x, h) = det Dh = (/in , ^Ta x ... x ^t„) 

where the cross product of spherical derivatives is controlled by using the 
point-wise Hadamard's inequality. 

(6.10) \hT,X ...xhTj^\hT,\---\hTj 
Hence 

(6.11) j(x,/i) ^ i/in| • i/iTr~^ 

Equality occurs if and only if the vectors h]^, hT2, hT„ are mutually or- 
thogonal, positively oriented and I/1T2I = ... = \hT„\- This amounts to saying 
that the Cauchy-Green tensor of h must be a diagonal matrix 

\hN\^ ••• 



(6.12) C{x,h)=D*h-Dh = 



Ihrf ■■■ 



••• |/iTp 



The notion of elasticity that we have introduced for radial stretchings can 
be extended to all weakly differentiable mappings. Let /i : A — )• A* be any 
mapping of annuli, and £ A be a point where h is differentiable.^ We say 
that, at this point, h is 

• expanding^ if | /lAr(xo) | > | hT{xo) \ , point of expansion 

• contracting, if | hN{xo) \ < \ hT{xo) | , Xo is a point of contraction 

In what follows we shall be concerned only with the regular points, that is 
II Dh{xo) II 7^ 0. We define the elasticity module of h at the regular points 
by the rule 

hN{x) I ^ ^ 



(6.13) 1] = r]h{x) = 



hrix) 



Homeomorphisms of class W^'"^ are differentiable at almost every point. 
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This is a measurable function taking values in the extended half line [0, oo]. 
Note the following inequalities, 

|D/t(a;)|" ^ |/ijvr 



(6.14) 



K.Q (x, h) 



J{x, h) 



I hN \ Ihrl 



n-l 



£)f^jn (^1 /lAT 1^ + (n- 1) I /tT 1^) 



n2j{x,f) 



n2 I h]\[ I I hr 



in-l 



(6.15) 



1 

Vh 



n-l 1 



n 



n 



We say that the deformation ^ : A ^ A* at the given point Xo is: 

• contracting, if rjh{xo) < 1 

• expanding, if rjh^Xo) > 1 

It will be clear later that the extremal deformations have the same elasticity 
type throughout the region. Furthermore, the elasticity module will not 
vanish inside the region. Thus for the extremal deformations the equation 
at (6.12) takes the form of a Beltrami system 

(6.16) D*h- Dh = J{x,h)nK{x) 

where the distortion tensor K is a diagonal matrix of determinant 1 



(6.17) 



K 



V 




















For the radial stretching 

(6.18) h{x) = 

we find that 

(6.19) hf^ = H{\x\)-N, 



Hi\x\) 



■Ti 



Hence, the cross product of the spherical derivative is parallel to the radial 
vector. 



(6.20) 



X ... X hT„ = { Tl ] n 



H 



n-l 



X 



\x\ 



The square of Hilbert-Schmidt norm of the differential is given by 

ij2 



(6.21) 



Dhf = + {n- 1) 



t = \x\ 



and the Jacobian determinant equals: 



n-l 



(6.22) 



J{x, h) 



HH 
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6.2. Some differential forms 

Integration of nonlinear differential expressions is best handled by using 
differential forms. Let us begin with the function t = \x\ and its differential 

(6.23) dt = ^idxi+ +XndXn 

\x\ 

The Hodge star of dt is an (n — l)-form. 

1 " 

(6.24) -kdt = -r—r '^{-Ifxi dxi A ... A dxi-i A dxj+i A ... A dxn 

' ' 1=1 

This form integrates naturally on spheres centered at the origin. We wish 
to normalize this form in such a way that the integrals will be independent 
of the sphere. In this way we come to what is known as standard area form 
on S"-i 

:*rdt xidxi A ... A dxi-i A dxi-^-i A ... A dxn 



(6.25) coix) = = j;(-l) 



i=l ' ' 



Viewing a; as a differential form on punctured space , we find that du = 0. 
By Stokes' theorem 

(6.26) / oj = ojn-i , for alH > 

J\x\=t 

Next, given any mapping h = {h^, h") iQ^Wlof class W^l'^{n,M.^), we 
consider the puUback of co via h 

(6 27) h^uj - \"(-lY^'^^^ ^ ■■■ ^ ^ ^^'^^ ^ ■■■ ^ ^^"^ 

1=1 ' ' 

Under suitable regularity hypothesis, for instance if /i G ■^oc" ^(^)I^")) 
and \h{x)\ ^ const > 0, this form is also closed, meaning that d{h*oj) = 
h* (duj) = 0. We view the exterior derivative of h as an n-tuple of 1-forms, 
dh = [dh^, d/i"] . With such a view dhAdt becomes an n-tuple of 2-forms. 
Further notation is self explanatory. For example, 

(6.28) \dh A dtp = \dh^ Adtf + ... + Idh"- A dt\'^ 
Now the decomposition at (6.8) takes the form 

(6.29) II Dhf = |d/ip = \dh A i^dtf + \dh A dt\'^ 

The reader may verify this formula by the following general algebraic identity 
|a|2 = |a A ★Afp + jo A iVp for ah covectors a € A^(l^'')- Just apply it to 
a* = d^*, i = 1,2, ...,n and N = dt, then add the resulting identities. A 
moment of reflection about (6.29) reveals that 

(6.30) dh A *dt = Hn dx 
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and 
(6.31) 



dh Adt = ^ {xihxj — Xjhxi) dxi A dxj 



Hence 

(6.32) \dh A dtf = [n - l)\hT\'^ 
Proposition 6.1. The following point-wise estimates hold 

(6.33) |/iAr| ^ A*di| 



(6.34) 



\hTr' > \h\ 



n-1 



dtAh^ 



CO 



Proof. Formula (6.30) holds for scalar functions in place of h we apply 

it to \h\ 



(6.35) 



|/iiv| ^ I \h\N I = I d\h\ A -kdt I 



establishing the inequality (6.33). 

The proof of the estimate (6.34) is more involved. Using (6.32), we are 
reduced to showing that 
(6.36) 

^{-ly—dh^ A ... A dh'-^ A dh'+^ A ... A dh"" A dt 
By Schwarz inequality it suffices to prove that 



dh A dt 


n-l 


n-1 







(6.37) 



dh A dt 



n — 1 



n-1 



Note that here in both sides only spherical derivatives of h are significant; the 
terms containing hjy dt vanish after wedging them with dt. This observation 
permits us to replace dh^ by the covectors Oj = dh^ — hj^dt. We view Oj as 
elements of the space /\^(M"~^). Once this interpretation is accepted, the 
proof continues via an algebraic inequality. 

Lemma 6.2. For everyn-tuple of covectors inW^~^, ai, On G /\^{R'^'~^), 
we have 



(6.38) |aip + ... + |a„p ^ (n-1) ( ^1 1^*1 - ^ ^i-^ ^ ^ - ^ ^\ 



1 

n-1 



Proof. We look at the matrix M of size n x (n — 1) whose rows are 
made of the covectors ai,...,a„. Let b2,...,bn € IR" denote the column 
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vectors of M 



(6.39) 



- - ai - - 

- - 02 - - 

- — Cln — — 



M 



bs • • • br 



In the left hand side of (6.38) we are actually dealing with the square of the 
Hilbert-Schmidt norm of M 



\ai\'^ + ... + |a„p = Tr (M*M) = jbsl^ + ... + |b„p 



(6.40) 

On the other hand, the sum in the right hand side of (6.38) is the Hilbert- 
Schmidt square of the norm of the cross product of the column vectors 

n 

^ |ai A ... A Oi-i A Oi+i A ... A a„|^ = |b2 x ... x bnf ^ Ihl^ • • • \K\^ 

i=l 

by Hadamard's inequality. The inequality between arithmetic and geometric 
means yields 

(6.41) " '2 " '2 - ' +- + K\ 



|b2|'---|bn|' ^ 



n- 1 



Prom here Lemma 6.2 is immediate. This also completes the proof of Propo- 
sition 6.1. 



CHAPTER 7 



Free Lagrangians 

Wc arc about to introduce the basic concept of this paper. Given two 
domains 0, and fi* in W^, we shah consider orientation preserving homeo- 
morphisms /t : $7 — >■ fi* in a suitable Sobolev class /^^'^(ri,^*) so that a 
given energy integral 

(7.1) I[h] = [ S{x,h,Dh)dx 

Jn 

is well defined. The term free Lagrangian pertains to a differential ra-form 
£{x, h, Dh) dx whose integral depends only on the homotopy class of h : 
$7 Vt*. The Jacobian determinant is pretty obvious example; we have, 

(7.2) I J{x,h)dx = \n*\ 

Jn 

This identity holds for all orientation preserving homcomorphisms of Sobolev 
class ri*). Many more differential expressions enjoy a property such 

as this. In the next three lemmas we collect examples of free Lagrangians 
for homeomorphisms h : A A* of annuli A = {x ; r < |x| < R} and 
A* = {x; r* < |x| < i?*}. 

Lemma 7.1. Let $ : [r*,i?*] — >■ R 6e any integrable function. Then the 
n-form 

(7.3) $(|/i|)d/i^ A...AcZ/i" 
is a free Lagrangian. Precisely, we have 

(7.4) / ^i\h\)J{x,h)dx = ujn-i I *T"-^$(r)dT 

for every orientation preserving homeomorphism h G W^''^{A,A*). 

This is none other than a general formula of integration by substitution. 
Lemma 7.2. The following differential n-form 

Xi dxi A ... A dxi-i A d\h\ A dx^+i A ... A dx„ (<^|^|) A -kdt 



(7-5) 



i=l 



\h\\x\^ \h\t 



n-l 



is a free Lagrangian in the class of all homeomorphisms h G W^'^(A,A*) 
that preserve the order of the boundary components of the annuli A and A* . 
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Precisely, we have 

Proof. First observe that the function : A — )• (r*,i?*) extends con- 
tinuously to the closure of A so that \h{x)\ = r* for |x| = r and \h{x)\ = 
for |x| = R. The point is that our arguments will actually work for arbitrary 
continuous mappings h : A ^ MJl of Sobolev class ^^'^(A, M") which sat- 
isfy the above boundary conditions, not necessarily homeomorphisms. The 
advantage in such a generality is that we can apply smooth approximation 
of h without worrying about injectivity. This extends Lemma 7.2 to all per- 
missible mappings. If h is smooth it is legitimate to apply Stokes' theorem. 
Recalling the form co and (6.26), we write the integral at (7.6) as: 



/ 

Ja 



{d\og\h\) Aoj = / d{ujlog\h\) = (logi?*) / a;-(logr*) / u 

JA J\x\=R J\x\=r 



(7.7) = con-i log — = Mod A* 
The same proof works for more general integrals of the form 

(7.8) ^ ( \h\ ) = u^n-i [m*) - Hn)] 
with $ G "^-"^[r*, i?*]. This also reads as 

J A Fl 

Another example of a free Lagrangian rests on the topological degree. 
Let * : 5""^ — >■ S"~^ be any smooth mapping of a sphere S'^~^ = {x ; \x\ = 
t} into the unit sphere C M". Then there is an integer, denoted by 

deg^*, and called the degree of such that 



(7.10) / *«(a;) 



a;„_ideg* 



where u stands for the standard surface measure of the unit sphere S"""^ C 
M", as given in (6.25). Note a general fact that deg^* remains unchanged 
under small perturbations of If ^ is an orientation-preserving diffeomor- 
phism, then dcg'I' = 1. Now a free Lagrangian somewhat dual to that in 
Lemma 7.2, is obtained as follows. 

Lemma 7.3. The following differential n-form 

, J^h' dh"" A ... A dh^-^ A d\x\ A dh^+'^ A ... A d/i" 

g \AW 
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is a free Lagrangian in the class of all orientation preserving homeomorphism 
mappings h G W^''^~^{A,A*). Precisely, we have 

(7.12) J J A h^uj = Mod A 

Proof. As a first step, we approximate h by smooth mappings h^ con- 
verging to h c-uniformly on A and in the norm of W^''^~^{A,W^). For v 
sufficiently large we consider the mappings 

(7.13) ^ = _^:A'^§"~\ u = l,2,... 

defined on a slightly smaller annulus A' = { a; ; r' ^ |a;| ^ i?'}, where 
r < r' < R' < R. The degree of ^ restricted to any concentric sphere in A' 
is equal to 1. Moreover, we have the point- wise identity '^'^oj = hi{oj). Now 
we integrate as follows 

(7.14) / —a¥uj=[ d(log\t\¥ 

Ja' i J A' ^ 

because d {^'^u) = 5'" (du) = 0. By Stokes' theorem, the right hand side is 
equal to: 

/ log |x| ^'^o; - / log |x| ^'''w = Wn-ilogi?' - Wn-ilogr-' = ModA' 

J\x\=R' J\x\=r' 

Passing to the limit as — > oo we conclude with the formula 

(7.15) [ — A^«a; = ModA' 

J A' t 

Finally letting r' ^ r and R' ^ R yield (7.12). 

Remark 7.4. The reader may wish to generalize this lemma to mappings 
/i : A ^ M" which are not necessarily homeomorphisms. The degree of J^ 
will emerge as a factor in front of the right hand side of (7.12). 

The following corollary from Lemmas 7.1, 7.2, 7.3 provide us with three 
inequalities, which we call free Lagrangian estimates. 

Corollary 7.5. Let h be a homeomorphism between spherical rings A 
and A* in the Sobolev class W^^'^{A,A*). Then 

(7.16) / ^\h\) \hN\ I/itI"-' > t^-^^{t) dr 

J A Jr, 

whenever is integrable in [r*,i?*]. We have the equality in (7.16) if and 
only if\hN\ I^tI""^ = J{x,h). Furthermore, 



(7.17) / 

Ja 

(7.18) / 

Ja 



\h\ \x 



|n-l 



> Mod A* 



m-i 



> Mod A 



l^l'*-! \x\ 
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Note that we have equahties if /i is a radial mapping. Other cases of 
equaUties are also possible. 



CHAPTER 8 



Some estimates of free Lagrangians 

Before we proceed to the general proofs, it is instructive to take on some 
estimates which are at the heart of our apphcation of free Lagrangians. 

8.1. The ^/i-energy integral with operator norm 

The extremal problem is remarkably simpler if we use the operator norm 
of the differential. Let hhe a permissible map of A onto A*. The following 
point-wise inequality is straightforward, 

(8.1) \Dh\ ^ max{|/i7v| , |^t|} 
Indeed, we have 

(8.2) \Dh\ ^ max{|/i7v|, I^Tal, IhrJ} ^ max{|/iAr| , {hrl} 

Certainly, equality occurs when /i is a radial stretching. All that is needed 
is to apply the lower bounds at (7.17) and (7.18). 

^ ' J A l^h I A l^h J A l^h 

Here we estimate each term by Holder's inequality. 



[Ml > (f 

Ja \h\^ ^ \Ja 



l^lkl""V \Ja\x\'', 
(8.4) ^ (ModA*)"(ModA)"-^ = a"ModA 

where 

def Mod A* 

^ ' ModA 
Similarly, we obtain 



1 

l-n 



[Ml > f f i^^i"'' ^ ( f 

Ja \h\^ " \JA\h\^-'\x\J UaNV 

n 1 

(8.6) ^ (Mod A) ~ (ModA) ~ = ModA 

Substituting these two estimates into (8.3) we conclude with the following 
Theorem 
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Theorem 8.1. Let A and A* be spherical rings in M", n ^ 2. Then for 
every h E ^{A,A*) we have 

(8.7) / ^ max{l,a"}ModA 

As for the sharpness of this estimate we note that equality holds for the 

power stretching h{x) = \x\"'~^x. 

There are, however, other cases of equahty in (8.7) if Mod A* ^ Mod A. 

8.2. Radial symmetry 

Suppose that the radial stretching 

(8.8) h{x) = H{\x\)^^, Heafe[r,R\ 

maps homeomorphically an annulus A = {x ; r < < R} onto A* = 
{y; < \y\ < i?*}, where ModA* = aModA, with some a > 0. We 
may assume that = and i?* = i?". Thus, in particular, h preserves 
the order of the boundary components. The power stretching h" = \x\"^^x 
serves as an example of such a homeomorphism. We shall prove that /i" is 
in fact the minimizer of among all radial stretchings. 

Proposition 8.2. For each radial stretching we have 

D/i ir . f w Dh"' ir . 2 



Equality holds modulo rotation of A. 

Proof. We express the energy of h as 

(\hj,\^ + {n-l)\hTf)^ 



i.W) ^h= f 

Ja 



\h\- 

where the radial and spherical derivatives of h are found in terms of the 
strain function H by using formulas (6.19); 

(8.11) |/tAr| = |i^(t)| and = 

Integration in polar coordinates leads to a line integral 



f^ft'^H^ \' dt 



R f 2. 

(8.12) = uJn-ilog-4 + dn 

where the integral average is taken with respect to the measure d/x = ^ 



The key is that X{t) = ^ is a one dimensional free Lagrangian. Indeed, we 
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see that the integral mean 



(8.13) 



-[ Xdfi 

Jr 



1 Hit) dt log I 



log I 



a 



does not depend on H. Next using Jensen's inequality for the convex func- 
tion X {X^ + ra - 1)5 yields 



+ n- 1 



Mod A 



(8.14) 



[a^ + n- 1] 2 Mod. 



as desired. Examining these arguments backwards we obtain the uniqueness 
statement. 

As a corollary, letting a go to 0, we observe that among all radial defor- 
mations h : A S"~^, the n-harminic energy 



(8.15) 



A= f 
JA 



\Dh\\'' ^ (n- 1)2 Mod. 



assumes its minimal value for the mapping h{x) = j^. This will be later 

shown to be true in a larger class of mappings. 

In dimensions n = 2,3 the assertion of Proposition 8.2 remains true for 
all permissible mappings in the Sobolev class #'^'"'(A, A*).^ 



8.3. Proof of Theorem 1.14 



The proof in case n = 2 is accomplished by an elegant use of Holder's 
inequality, 



(8.16) 



Dh\ 



\h\' 



flN 



{Llm) +(X 



+ 

hr I 



\x\ \h\ 



dx 



\x\ 



The lower bounds at (7.17) and (7.18) yield 



as claimed. 



(ModA*)^ + (ModA)2j (ModA)"^ = (q^ + l) ModA = 



In the next chapter we shall extend this result to all dimensions, but it will require 
an upper bound of the modulus of A* . 
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The proof in dimension n = 3 involves application of Holder's inequality 
twice 



J 

JA 



Dh\ 







2 

+ 2 




• 




h 


h 






dx 



Now using lower bounds at (7.17) and (7.18) we find that 



(ModA*)^(ModA)~^ + 2ModA ' (ModA)-t 



+ 2)2 Mod A 



(8.17) 

For the uniqueness we refer to Section 9.3 where such problems are dealt 
in all dimensions. The lower bounds of free Lagrangians, for the energy 
integrand j^^J in higher dimensions, are more sophisticated. We found 
optimal bounds if Mod A* ^ a^ModA, where 1 < a„ < oo is a solution to 
the algebraic equation 



(8.18) 



(a^ + n - 1)" (a^ " l) = "n" ^ < 



n 



n ■ 



CHAPTER 9 



Proof of Theorem 1.15 



It suffices to consider homeomorphisms /i : A — t- A* of class #'^'"(A, A*). 
Recall that the target anniihis is not too fat. Precisely we assume that 

(9.1) ModA* < a„ModA 

where the critical factor a„ > 1 is determined by the equation (8.18), which 
we write as 

Lemma 9.1. Let X,Y ^0 and 1 ^ a < a„. Then 

Q^ f \ def (g^ + 71-1)^^(0^-1) 

(9.3) a = a(a) = < 1 

and, we have 

(9.4) [X^ + (ra - l)y2] ^ ^ a X" + 6 XY""^ 
where 

n-2 

n [o? + n — l) 2 

9.5) h = ^ 

a 

Equality holds if and only if X = aY. 

Proof. Because of homogeneity we may assume that Y = 1; the case 
y = is obvious. We are reduced to proving an inequality with one variable 
X^O, 

n 

(9.6) [X^ + (n - 1)] 2 ^ aX" + bX 
Consider the function 

(9.7) (p{X) = (X^ + n - 1) 2 - aX"" - bX 
and its two derivatives 

(9.8) ip'{X) = n {X"^ + n - l)"^ X - naX""-^ - b 



(9.9) ip"{X) = n(n-l) 



n — 4 



{X^ + n-1) 2 {X^ + l)-aX'' 



> 



This letter estimate is guaranteed by our hypothesis that a < 1. Thus ip is 
strictly convex and, therefore, has at most one critical point. The coefficient 
b has been defined above exactly in a way to ensure that ip'{a) = 0. The 
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coefficient a has been chosen so that (p{a) = 0. This completes the proof of 
Lemma 9.1. 

Lemma 9.2. Let X,Y ^ and ^ a ^ 1. Then 

(9.10) [x2 + (n-l)y^]* ^aF" + 6Xy"-^ 
where 

(9.11) a = (n- 1) (a^ + n- l) ^ (l - a^) 

and 

(9.12) b = na {a"^ + n - l)"^ 

Equality holds if and only if X = aY. 

Proof. Because of homogeneity we may assume that Y = 1, the case 
y = is obvious. We are reduced to the inequahty with one variable X ^ 0, 

n 

(9.13) [X^ + (n - 1)] 2 ^a + bX 
Consider the function 

(9.14) ip{X) = {X^ + n - l)'^ - bX 
and its two derivatives 

n-2 

(9.15) (p'{X) = n{X^ + n-l) ^ X-b 

n—4 

(9.16) ip"{X) = n{n-l){X^ + n-l) ^ (X^ + 1) > 

As before, (p is strictly convex with one critical point. The coefficient b has 
been defined exactly to ensure that (^'(a) = 0, whereas a = ip{a). This 
completes the proof of Lemma 9.2. 

We now proceed to the proof of Theorem 1.15. Consider two cases. 

9.1. The case of expanding pair 

We apply Lemma 9.1 with X = y = ^ and a = > 1, to 

obtain the point-wise inequality 



+ (n - 1) 



(9.17) ^ o— r; \- b m 

Then we integrate (9.17) over the annulus A. For the last term we may 
apply the lower bound at (7.16). The first term in the right hand side of 
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(9.17) needs some adjustments before using the lower bound at (7.17). These 
adjustments are easily accomplished by Holder's inequality. 



/ 

JA 



\h\^ \JA\h\\x\^-^J \Ja\x\^J 



^ a(ModA*)"(ModA)^-" + 6ModA* 
= (aa" + 6a)ModA 

(9.18) = (a^ + n-1)^ ModA 

because we have equality at (9.6) for X = a. We see that the right hand 
side is the energy of the power stretching 

(9.19) /i"(x) = rV-"|x|"-^x, a 



Mod A 



9.2. The case of contracting pair 

This time we use Lemma 9.2 with the same data X, Y, a as were used 
in the case of expanding pair; that is, with X = -l^, Y = ^-j^ and a = 

^^J^ ^ 1. This gives the following point wise inequality 
(9.20) ^ + ^'-'^^' 




We integrate it over the annulus A. For the last term we apply the lower 
bound at (7.16). As in previous case the first term in the right hand side 
of (9.20) needs some adjustments before using the lower bound at (7.18). 
These adjustments are easily accomplished by Holder's inequality. 



^ a (Mod A) ^ (Mod A)^ +b Mod A* 
= (0 + 60;) Mod A 

n 

(9.21) = (a^ + n- 1)2 Mod A 

because we have equality in (9.13) for X = a. We easily see that the right 
hand side is the energy of the power stretching 

(9.22) /i"(x) = rV-"|x|"-^x, a 



Mod A 
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9.3. Uniqueness 



The reader might want to compare this proof with Chapter 13, in which 
the uniqueness question is dealt in greater generahty. Let /i : A — >■ A* be any 
extremal mapping. First observe that in both cases we have used a general 
inequality 

(9.23) J{x,h) ^ |/ijv| Ihr^'^ 

For h to be extremal, we must have equality. In view of (6.11) the vectors 
hN,hT2, ■■■,hTn are mutually orthogonal and I/itjI = ••• = |/it„| = Ihrl- 
Using the matrix representation of Dh at (6.5) we find the Cauchy-Green 
tensor of /i to be a diagonal matrix 

|2 



(9.24) 



D*hDh = 



\hN\' 














\hj 



Another necessary condition for h to be extremal is that |/iAr| = a|/ir|, 
because of the equality cases in Lemmas 9.1 and 9.2. In this way, we arrive 
at the Beltrami type system for the extremal mapping 



(9.25) D*hDh = J{x,hy^ 

where K is a constant diagonal matrix 

a^-n 



K 



(9.26) 



K 







a 








a 





The power mapping /i" = h"'{x) is one of the homeomorphic solutions to 
this system. It is well known and easy to verify that other solutions are 
obtained by composing this particular one with a conformal transformation 
[7, 43, 44]. Thus h takes the form 

(9.27) h = goh'^ 

where g' : A* — )• A* is a conformal automorphism of the target annulus 
A* = 74(r*,M*) onto itself. Thus, up to rotation ^(y) = y oy g{y) = r^R^-^. 
In either case /i is a power stretching up to the rotation. The proof of 
uniqueness is complete. 



Part 2 



The n-Harmonic Energy 



A study of the extremal problems for the conformal energy 

(9.28) ^h= [ \\Dh{x) W'^dx 

Ja 

is far more involved. Even when the solutions turn out to be radial stretch- 
ings, they no longer represent elementary functions such as = \x\"~^x. 
And this is not the only difficulty; there are new phenomena in case of 
non-injcctivc solutions. In view of these concerns, it is remarkable that the 
method of free Lagrangians is still effective for non-injective solutions. First 
we examine planar mappings. 

Concerning mappings minimizing the ^^-norm of the gradient, we ref- 
eree the interested reader to [18] and references there. 



CHAPTER 10 



Harmonic mappings between planar annuli, 
proof of Theorem 1.8 

Let /i : A — ^ A* be a homeomorphism between annuli in the Sobolev 
class W^''^{A, A*). We view h as a complex valued function. Let us recall 
the formulas 

(10.1) lDh\\^ = \hNf + \hT\^ 

and 

(10.2) det Dh = Im {hrh^) ^ |/iiv| {hrl 

Case 1. The expanding pair, ^ -pr- 

Thus the target annulus is conformally fatter than the domain. We find 
a unique number w ^ such that 



(10.3) 5 = Sl±V^ 

r n + - Lo 

Without loss of generality we may assume, by rescaling the annulus A if 
necessary, that 

(10.4) R = + \/Rl-oj and r = r* + \/rl-u 
We begin with the inequality 

(10.5) f^iJM -imV^o 



Equivalently, 



Here we have used an elementary fact that | /iat | ^ | |/i|iv |, equality occurs 
if and only if ^ is a real valued function. Let us integrate this estimate. 
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We apply Holder's inequality to the first term and the estimate (7.16) to 
the second term, 

J \\Dh\\^ ^ ~^ j [log(l^l + vW^)^ 

f-R, ^2 

+ / 

Jr* V — W 



N \x\ 



dx 
A 



+ 



-U) 



27r log 



R* + s/ R* 



UJ 



, R 

ZTT log — 

r 



-1 



+ 



+ 27r 



— w + a;log(r + \J — u:) 



:io.6) 



OJ 



rf — w 



Elementary inspection reveals that equality holds for the Nitsche mapping 



1 



h{z) = -(z + - 



OJ 



Case 2. The contracting pair, < ^ ^ + 



R ^ R, 



R'i 



1. 



In this case the target annulus is conformally thinner than the domain, 
but not too thin. We express this condition by using the equation (10.3), 
where this time < a; ^ r^. Again we may assume that the relations at 
(10.4) hold. The same inequality (10.5) can be expressed in a somewhat 
different form 

2 



(10.7) 



DhV >LJ 



h 



As before, we apply Holder's inequality, which together with the estimate 
(7.16) yields 

2 



Dhr ^oj 



\x\ \h\ 

Next, the estimate (7.18) gives 



I 

J A 



Dh\ 



^ a; ( 27r log 



R 



27r log 



R 



-1 I'R, 

+ Att I \Jt — ijj dr 



(10.8) 



, R 
= ZTTLO log h 27r 



■ U — LJ 



l0g(T + \Jt'^ —Uj) 



In view of (10.4) we find that 



(10.9) 



/ \Dh\'^ ^:^2^lR^^/^l^ - 2'Kr^^/^ 

J A 



rt — u 



Again equality holds for the Nitsche mapping 

M.) 4(^ + 1) 



10. HARMONIC MAPPINGS BETWEEN PLANAR ANNULI 



73 



The borderline case. Taking w = r J we obtain what is called the critical 
Nitsche map with R = + \/R'^ — r1 and r = r* , 

(10.10) h\z) = \[z+'i^ 




Figure 1. The Jacobian of h* vanishes on the inner circle. 

The confer mal energy of the critical Nitsche map equals: 
(10.11) / II Dh* p = 

J A 

We arc now in a position to consider the case: 

Case 3. Below the lower Nitsche bound, ^ + -1 < ^. 

The target annulus A* is too thin. We shall see that an inner part of 
A has to be hammered flat to the inner circle of A*. We can certainly 
rescale the annulus A to have R = R^ + \/R'^ — r^. This together with the 
hypothesis of this case, yields r < . 



A 




The critical Nitsche To be squeezed onto 
annulus the inner circle 



Figure 2. Below the lower Nitsche bound. 



For every permissible map /t : A ^ A* we still have the estimate (10.8), 
in which we let oj = r'i. Hence 



Dh\ 



R 

^ 27rr* log h 27r 



TV — — r. 



log(r + Vr^ 



(10.12) 



= 2ttR, ^JRI - rl + 2T,rl log ■ 
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The first term represents conformal energy of h* : A(r*,i?) — )■ A{r^:,R^:), 
by (10.11). The second term turns out to be exactly the energy of the 
hammering mapping g{z) = J^* |f| , which takes the remaining part A(r, r*) C 
A onto the inner circle of A(r*, i?*). 



is smaller than the energy of any homeomorphism h: A(r, R) A(r*, i?*). 
It is clear that h° is a ^^'^-limit of such homeomorphisms, completing the 
proof of Theorem 1.8. 

Now we proceed to higher dimensions. Different estimates will be re- 
quired for the contracting pairs of annuli than for the expanding pairs. Thus 
we devote a separate section for each case. 




Thus the energy of the mapping h° : A A* 





CHAPTER 11 



Contracting pair, Mod A* ^ Mod A 



First we consider the case when A* is not too thin. 

11.1. Proof of Theorem 1.12 

Here we assume the Nitsche bound, A/f (Mod A) ^ Mod A*, see (1.35) for 
A/f. This bound means precisely that there is a radial n-harmonic homeo- 
morphism 

X 



(11.1) 



/i° : A 



h°{x) = H{\x\) 



\x\ 



Recall the characterictic equation at (5.28) for H = H{t) 



(11.2) 



n 



where c is a positive constant determined uniquely by the rings A and A* 
Equivalently, 



(11.3) 



1 + 



TC-2 

2 



n 



In particular, c ^ [if(r)]" = r". Equation (11.3) suggests that we should 
consider the nonnegative solution r) = r){t) to the equation 



(11.4) 



1 + 



V 



n — 2 
2 



^ 1 , r* < t < i?* 



n-lj ' ' ' t"" 
There is exactly one such solution. The values of rj lie in the interval [0, 1] . 

Now, let /i : A A* be any homeomorphism of annuH preserving 
both orientation and order of the boundary components. We assume that 
h G ^^'"(A,A*). For each a; G A, we apply Lemma 9.2 with X = \h^{x)\, 
Y = \hj.{x)\ and a = ri{jh{x)\) ^ 1, to obtain the point-wise inequality 



.Dh\ 



|/i,v|' + (n-l)|/ir|' 



(11.5) 



\ — 

^ (n — 1) 2 c 



M 

\h\^ 



+ b{\h\) \hN\ \hT\ 



n-l 



The coefficient b{\h\) comes from (9.12) where we take a = ?7(|/i|). An 
important fact about b = b(jh\) is that we have equality at (11.5) if |/iAr| = 
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r/(|/i|) |/it|- This is exactly happening for the radial n-harmonic map at 
(11.1), by the definition of the constant c. Let us integrate (11.5) over the 
ring A. For the last term we may apply the lower bound at (7.16). For the 
first term in the right hand side of (11.5) we use Holder's inequality, and 
then (7.18). 



/ 

Ja 



(11.6) +con-i I *T"-^6(r)dT 

rti, 

^ (n- l)icModA + Un-i I T"-~'^b{T) 

Jr» 



dT 



Finally, observe that we have equalities in all estimates for the radial stretch- 
ings. Thus 



(11.7) [ II ^ / 

J A J A 



Dh° 



as claimed. A proof of the uniqueness statement is postponed until Chap- 
ter 13. We only record further use that the equality at (11.7) yield 

(11.8) |/ijv| I/itI""^ = J{x,h) 

Next we turn to the case when A* is significantly thinner than A. 

11.2. Proof of Theorem 1.13 

Here we assume that the lower Nitsche bound fails; that is. 

Mod A* <A/'t(ModA) 
We split A into two concentric annuli 

(11.9) A = A(r, i?) = A(r, 1) U A[l, i?) 
where 

(11.10) ModA* = A/'t(ModA[l,i2)) 

Let H : A[l,i2) — )• A* denote by critical Nitsche map. This is the radial 
n-harmonic function 



(11.11) ^x) = H{\x\)-^^ 

determined by the characteristic equation 



(11.12) 



2 t'H' 



n 
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Evaluating it at t = 1 yields c = = . Now consider any permissible 

map /t : A ^ A* in ^(A, A*). We may use inequality (11.6) with c = to 
obtain 

(11.13) / ||m|r >(n- 1) ir^ Mod A + a;„_i /" *T"-^6(r)dT 

J A Jrt 

On the other hand we have the following identity for the critical Nitsche 
map 

(11.14) / ||DH||" = (n-l)trJ^ModA[l,i2)+w„_i /^*T"-i6(T)dr 

Ja[1,R] Jr» 

Hence (11.13) takes the form 

(11.15) f ||m|r ^ (n-l)tr:^ModA(r,l)+ /" || -Di^ 

As a final step we notice that the first term is precisely equal to the energy 
of the hammering map 

(11.16) (n- l)5r^Mod A(r,l) = / \\Dg\\'' , g{x) = r^-^ 

A(r,l) FI 

We now glue K and g along the sphere |x| = 1 to obtain a map h° : A ^ A* 

(11.17) h° 



] g{x) on A(r, 1) 
\)^{x) onA[l,ii:) 

This map minimizes the conformal energy. Indeed, 
(11.18) / \\Dhr^ [ ||m°||" 

J A J A 

for every h G ^{A,A*). The question of uniqueness is discussed in Chapter 
13. 



CHAPTER 12 



Expanding pair, Mod A* > Mod A 



Such annuh determine uniquely a radial n-harmonic map h° : A ^ A* 

in the class H^, 

(12.1) h°{x) = H{\x\) where CH = c<0 

If A* is too fat, it will latter become clear that in dimensions n ^ 4 this 
radial n-harmonic mapping is not the minimum energy solution. 

12.1. Within the bounds, Mod A < Mod A* < A/'^(Mod A) 

Here the upper Nitsche function J\f^ = J\f\t), t > 0, is determined 
uniquely by requiring that the inequality Mod A* ^ J\f^ (Mod A) be equiva- 
lent to 



■n — 

(n'^ + n-l) 2 („2 _ i\ 

(12.2) ^ ^-^ ^ ^ 1 

where rj^ = ^ is the elasticity function of : A — t- A*. A somewhat 
explicit formula for is given in (12.13) after an analysis of (12.2). It may 
be worth noting in advance that J\f^ = oo in dimensions n = 2,3. Thus 
the condition Mod A* ^ A^t(ModA) is actually void for n = 2,3. By the 
definition of Q!„ at (9.2), condition (12.2) is equivalent to 

(12.3) r]^ (t) ^ a„ , for all r ^ i ^ i? 

Since we arc in the expanding case the function rj^j is decreasing. We there- 
fore need only assume that 

(12.4) l<r]„{r)^an 

Recall that a2 = as = oo, so this condition poses no restriction on A* in 
dimensions n = 2 and 3. 

Although it is not immediately clear the condition (12.4) is in fact a 
condition on the moduli of A and A* alone, precisely 

(12.5) F_(7n)-^^i?-(7nf) ' 



^In dimensions n = 2,3 we have a„ = ca, hence 7n = 1. In this case H-{'y„) = 
and, therefore, we impose no upper bound for Mod A*. 
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(3^+„-i)^(3^-i) 



Figure 1. The bounds for the elasticity of h°. 

where the number 7„ > 1 (for n ^ 4) is determined by the equation 
(12-6) 7n = r_ (^) 

Proof of (12.5). Since h° hes in the class %- its strain function takes the 
form 

(12.7) H{t) = XH(kt), k>-, for r^t^R 

r 

The boundary constrains H(r) = and H(R) = yield a system of 
equations for the parameters A and k 

(12.8) \H_{kr) = r^ and XH^{kR) = R^ 

Once we eliminate A the parameter k is determined (uniquely) from the 
equation 

(129) R._ium 

To solve this equation we look at the function 

(12.10) Q{k) = , > 1 
see (5.64). Elementary computation shows that 

(12.11) =^_(i^A;)-r?_(rfc)<0 

In view of (5.65) we see that Q{k) decreases from +oo to ^ as runs from 
^ to GO. Recall that we are in the expanding case, > 7 ■ Then by Mean 
Value Theorem there is unique k > ^ satisfying (12.9). We now rewrite 
condition (12.6) as follows: 

77^ (r) ^ an <^ r]_ {k r) ^ <^ 

(12.12) kr ^ n:\an) = r^{^)=jn ^ 
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Here we have used the identity r_(s) = 7!^ (i), see (5.49) and (5.53). The 
latter inequaUty is equivalent to Q (^) ^ Qik) = j^, which is the same as 
(12.5). We are now in a position to define 

(12.13) M\t) = UJn-l logH_ ( 7„exp ^ ) - Un-l log//-(7n)- 

Proposition 12.1. Let h -. A ^ A* be a permissible map in ^(A,A*) 
where 

(12.14) 5 < ^1 < (~'»f ) 



Then 



(12.15) / \\Dhr^ [ 

J A J A 



Dh° 



Equality holds if and only if h{x) = h°{x) up to a conformal automorphism 
of A. 

Proof. The characteristic equation (5.28) for the map h° : A ^ A* 
defines a positive constant q = —{n — 1)~2~ JCH; that is, 

(12.16) q = t^H^ + (n - 1)H^ ' (t^H'^ - H^^ > 
Equivalently, 

n-2 

(12.17) ^n-l + vl) ^ {r,l-l) = ^^ 



Now, consider an arbitrary permissible map /i : A — > A*. We may assume 
that h is a homeomorphism preserving orientation and the order of the 
boundary components of the annuli. We introduce a function rj = r]{t) 
implicitly defined for r* ^ i ^ i?* by the equation, 

(12.18) („_l+^2)^(^2_-L)^^ 

Note that the equation (12.17) for rjjj &t t = r coincides with (12.18) for 
at t = r*. Hence ry(r*) = r]^{r). Since rj = ri{t) is decreasing we find that 
for every r* ^t^R^, 

(12.19) l<r,{t)^ri{r^) = r,„{r)^an 

Next we recall the function a = a(a), 1 ^ a < oo, from the formula (9.3). 
By the definition of a„ at (8.18) it follows that 

(12.20) a{r]) = (^"l + ^^J^^ (^^ - 1) ^ r^ < t < 



82 
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A 



a{\h\) 



1 '7(l''l) a„ 



Figure 2. The coefficient a{\h\). 



Now Lemma 9.1 applies to X = \hN\, Y = \hT\ and a = r]{\h\), 

n 



Dh\ 



(12.21) 



^\hNr + in-l)\hT\ 
^ a{\v\) l/iivr + KH) \hN\ \hT\ 
\hN\ 



n—l 



\h\vi\h\) 



+ ^{\h\)\hN\ \hT\ 



n-1 



wfiere = b(r]{\h\)^ and b is given by (9.5). According to Lemma 9.1, 

equality holds at a given point x if and only if I ft. AT (x) I = ?7(|/i(x)|) \hTix)\. In 
particular, it holds almost everywhere for h = h°{x), because \ h°j^\ = ry// 
We now integrate over the annulus A. The last term at (12.21) is easily 
handled by the lower bound at (7.16), 



'^{\h\) \hN\ \h 



in— 1 



(12.22) 



f *T"-i$(T)dr 

J r* 

= [ mh°\) \h%\\h°rr' 

J A 



The first term in the right hand side of (12.21) will be treated by Holder's 
inequality in order to apply the lower bound at (7.17). This should be done 
in a way so one obtains equality for h° . Taking into an account the identities 



(12.23) 

We proceed as follows 

\hN\ 



N\ 



H 



1 



\h°\ rj{\h°\) Hri^ \x\ 



(12.24) 



/ 

J A 



\h\ V{\h\) 



A \h\ ri{\h\) \x 



Ua\xy) 



l-n 



The first term in the right hand side is none other than a free Lagrangian 
at (7.9). This allows us to quickly evaluate the first term, without getting 
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buried under massive computation. 

\h\j^ dx 



r \h\j^ dx _ r \h°\,^dx _ r H{\x\)dx _ r 
J A \h\ 7j(\h\) " J A \h°\ v(\h°\) " A Hri^\x\^-^ ~ J a 



dx 



Ia \h\ 77(|/i.|) 
Hence 

(12.25) 

In conclusion, 



\hN\ 



\h\ r,{\h\) 



dx 



Mod A 



(12.26) / II D/i ir ^ g Mod A + a;n-i / *r"-^$(r)dT 

J A Jr» 
with equality occuring for h° , as claimed. 

The uniqueness results will be proven in more unified manner in Chap- 
ter 13. 

12.1.1. Proof of Theorem 1.10. Indeed, the condition Mod A* ^ 
A/'^(Mod A) is equivalent to (12.14), in view of the explicit formula for 
given in (12.13). 

12.1.2. Proof of Theorem 1.9. The proof is immediate from Propo- 
sition 12.1 once we recall that N"^ = oo, for n = 2, 3 posing no upper bound 
for Mod A*. 



CHAPTER 13 



The Uniqueness 



There is a comprehensive approach to all our questions concerning unique- 
ness of the conformal energy 

(13.1) <^h= I \Dh{x)\'^dx -minimal mappings /i : A ^ A* . 



We have already found that within the Nitsche bounds 

(13.2) A/"t(Mod A) ^ Mod A* ^ M^ModK) 
the radial n-harmonics 

(13.3) h°{x) = H{\x\) ^ 

are among the extremal solutions. We shall show that 

Theorem 13.1. Below the upper Nitsche bound; that is, 

(13.4) Mod A* ^ Af^ (Mod A) 

every permissible minimizer h : A A* coincides with the radial extremal 
map modulo conformal automorphisms of A. 

Let h G =^(A, A*) be any permissible extremal mapping. In all the 
preceding cases we came to the following equation as one of the necessary 
conditions for h to minimize the energy 

(13.5) \hM\=v{\h\) \hT\ 

see (9.17), (9.20), (11.5) and (12.21). Here ry = r]{T) is anonnegative function 
defined for r* ^ r ^ i?* by the rule 

(13.6) r'^ (77^ + n - 1) (r/^ - 1) = c 

Here the constant c comes from the characteristic equation for the radial 
extremal map h°{x) = H(jx\) Precisely, we have 

(13.7) c ^ [Hitr ivl+n-l)^ {nl - 1) , where r,, (t) = 

Note that rj is strictly positive if c > 0. However, if c < then r]{T) is 
strictly increasing and, therefore, can vanish only at the endpoint t = r^f. 
This latter situation arises when t]^ = in the hammering part of h°. 
Another necessary condition for a mapping to be extremal takes the form 

(13.8) |/tjv| l/iTl""^ = Jix,h) 
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because in each case we used the estimate (7.16) , which is sharp only 
when (13.8) holds. Consequently, the Cauchy-Green tensor of h becomes a 
diagonal matrix, 

l2 



(13.9) 



D*h-Dh = 



















\hT\' 



Mc(x) 



see (6.12). Whenever Holder's inequalities were used we always arrived at 
one of the following conditions for h to become extremal 



(13.10) 



hr I = —TT^ , or I /ijv I = — r-j-^r?(|/i|) 



\x\ 



where Ci and C2 are constants, see (11.6) and (12.24). No matter which 
case we look at, the conclusion is that 



(13.11) 



X\h\ 



and I hisf 



\\h\ 



v{\h\) 



for some constant A > 0, because of (13.5). This constant A is the same 
for all extremal solutions; in fact, A = 1 . This can be observed as follows. 
Since equality holds in either (11.6) or (12.25), at least one of the integrals 



/ 

Ja 



m-i 



\x\ \h\ 



n-l 



and 



/ 

Ja 



\hN\ 



\x\ri{\h\) 



is equal to Mod A. In view of (13.11) this yields A = 1. Hence 

A2|^|2 



(13.12) 



D*h ■ Dh 



\x\ 



V'{\h\) 
1 
















Gix,\h\) 



We emphasize that the metric tensor G = G(x,r), viewed as a function of 
two variables x E A and r G (r*,ii!*), is the same for all extremals. A key 
step in establishing uniqueness is: 

13.1. The Point-Cauchy Problem 

We shall consider the Cauchy-Green equation 
(13.13) D*h{x) Dh{x) = G(x, \h\) 

for mappings h : Q, ^ K* defined in a domain C M" and valued in an 
annulus A* = A (r*, i?*) C M". Here the function 



(13.14) G:nx{r^,R^) 



{the space of symmetric 
positive definite matrices 
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is assumed to be ^°°-smooth. As for the regularity of the solutions, we 
initially assume that h is only continuous. It then follows from the formula 

(13.15) \\Dh\\^ = TTG{x,\h\) G ^{^^in) 

that h is actually locally Lipschitz continuous. One can look at O as a Rie- 
mannian manifold equipped with the positive definite element of arclength 

ds^ = Gij{x) dec* (g) dx^ where [Gij{x)] = G{x, \h{x)\) 

In this way /i : — )• A* becomes a local isometry with respect to this metric 
tensor on ft. At this point we may appeal to the well known regularity result 
due to Calabi and Hartman [8]. It tells us that if D*h-Dh G '^'''°'{n,Rf') 
for some < a < 1 and = 0,1,..., then h G ^''+i'"(17, M"). This result 
can be apphed repeatedly to infer that in fact h G ^°°(fi,A*). One more 
observation is that the system (13.13) is invariant under linear isometrics of 
the target annulus A* . Precisely, if h solves this system so does the mapping 
Th for every linear isometry T : M" M". A priori the system (13.13) 
may not admit any local solution. But if it does, the Riemann curvature 
tensor of G{x, \ h{x)\) must vanish. We shall take advantage of the classical 
computation of curvature. We express the second derivatives of h in terms 
of its first order derivatives 

Here F^^ are the Christoffel symbols; explicitly, 

/-,o17^ _ 1 v^^ii^ f ^Gjj dGjk dGjk 

(13.17) + 

1=1 ^ ■' 

where G^" are the entries of the inverse matrix to G, see for instance [29, 
p. 37]. It should be noted that the partials are computed in accordance 
with the chain rule; so in fact we have 

(13.18) |^ = G.,+G.|^ 

dxk dxk 

where Gx^. and Gr stand for the partial derivatives of the function (x,t) 
G(x, r). Thus, in particular, the Christoffel symbols depend linearly on Dh. 
Let us state the general form of the second order equations obtained in this 
way 

(13.19) D'^h{x) = $(x, \h\, Dh) 

where $ : O x (n,i?*) x M"^" — ^ j^nxnxn given "^"^-smooth function. 
This is in fact a quadratic polynomial with respect to Dh. It should be 
noted that (13.19) is frame indifferent; that is, it holds for Th if it holds 
for h . With these equations at hand we can now prove the following 
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Lemma 13.2. (Uniqueness in the point Cauchy problem) Suppose 
we are given two solutions h° and h to the Cauchy-Green equation (13.13) 
such that \h{a)\ = |/i°(a)| for some point a G fi. Then there is an isometry 
r : M" ^ M" such that 

(13.20) h{x) = Th°{x) for all X en 
Proof. It follows from the equation (13.13) that 

(13.21) D*h{a) Dh{a) = D*h°{a) Dh°{a) 
Therefore, there is an isometry T : — >■ such that 

Dh{a) = ToDh° (a) = D {Th°) (a) 

With a suitable choice of the isometry we may assume without loss of gen- 
erality that the first order derivatives of h and h° also coincide at a; namely, 

(13.22) \h{a)\ = \h°{a)\ and Dh{a) = Dh°{a) 

We are going to show that these two equations hold in the entire domain 
ri. Obviously, the set where (13.22) hold is relatively closed in Q,. Thus, we 
need only show that this set is also open. To this end we consider a small 
ball 

]B = {xGR"; \x-a\^e} 

center at a and contained in Q. We estimate the supremum norms of Dh — 
Dh° over B with e approaching zero. 

II I^h — Dh° II _^oo^B-j ^ e II D^h — D^h° \\ _5foo(B) 

= e II ^x, \h\, Dh) - $(x, \h°\, Dh°) II ^oo(B) 
=^ e II \h\ - \h°\ II ifoo(B) + e || Dh — Dh° \\ ^oo(b) 

Here and subsequently, the symbol =<; stands for the inequality with a con- 
stant independent of e. This constant varies from line to line. Noting that 

(13.23) e II \h\ - \h°\ ||^cx=(b) =^ e || Dh - Dh° ||^oo(b) , 

we conclude that Dh = Dh° and \h\ = \h°\ in B, as desired. 

The equation Dh{x) = Dh°(x) implies that h°{x) = h{x) — 2v, where 
V is a constant vector. This combined with the condition |/i°(a;)| = \h{x)\ 
yields that h{x) — v is orthogonal to v. If v was not zero the image of J7 
under h would be an (n — l)-hyperplane and, consequently, J{x, h) = 0. 
This is impossible for any solution to the equation (13.13), because 

(13.24) J(x, h) = VdetG(x, \h\) ^ 
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13.2. Proof of Theorem 13.1 

We may assume that the extremal map /i : A — )• A* preserves both the 
orientation and the order of the boundary components of the annuh. For, if 
not, wc compose h with a suitable conformal automorphism of A. By virtue 
of Theorem 1.7 any extremal map h is monotone and 

(13.25) A (r,, i?,) = A* C h{A) C A^ = A [r*, i?,] 

Thus, in particular, the set h^^{A*) is a connected subset of A. Wc also con- 
sider the annulus A(p, R) = {h°)-\A*) C A. In case Mod A* ^ 7Vt(Mod A) 
the annulus A(p,R) is the entire annulus A = A{r,R). However, when 
Mod A* < A/t(Mod A) the inner radius of the preimage (/i°)~^(A*) = A(p, R) 
is determined in such a way that 

|^°(x)| = r* if r < |x| ^ /9 

and 

r* < \h°{x)\ < R* a p<\x\ < R 

We shall soon see that h-^{A*) = A{p,R) for every extremal map h. To 
this effect let us consider the union 

n = h-^{A*)U A{p,R) C A 

which is connected because /i"^(A*) n A(p, R) / 0. This latter set actually 
contains all points near the outer boundary component of A. Precisely we 

have lim \h{x)\ = i?*. We now make use of the free Lagrangian identities 

\x\^R 

(7.4) 

(13.26) / %^ = Mod A* = / 

It follows from (13.12) that 

77(|/i(x)|) dx f r](\h°{x)\^ dx 



Ja W Ja 



\X\ 



The point to make here is that both Jacobians J{x,h°) and J{x,h) vanish 
almost everywhere outside Q. Now, by mean value property for integrals 
there exists a point a G such that 

(13.28) v{\h{a)\)=r,{\h°{a)\) 

Let U C $1 denote the set of all points a € Q for which (13.28) holds; U 
is certainly relatively closed. This set U is also open. To see this we first 
observe that 'rj{\h{a)\) > for all a G U. Indeed, if a G A{p,R) then 
|/i°(a)| > r*, whereas for a G h~^{A*) we have \h{a)\ > r^. It remains to 
recall that ?7(r) > 0, whenever r > r^,. In other words, the Cauchy-Green 
tensor of h is positive definite near every point of U. Now it is legitimate to 
appeal to Lemma 13.2. Accordingly, there exists a linear isometry T : ^ 
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M" so that h{x) = Th°{x) near the point a; thus, \h{x)\ = \h°{x)\ near this 
point. This shows that U is the entire domain We then have 

\h{x)\ = \h°{x)\>r^ for all a; GO 

Again by Lemma 13.2 we infer that, upon suitable adjustment via isometry 
of the target annulus it holds 

h{x) = h°{x) in 

As a matter of fact we have 

Indeed 

X G /i"^(A*) <^ h{x) G A* <^ h°{x) G A* <^ x G k{p,R) 

Thus the uniqueness theorem is proven if A(p, R) = A. This is the case 
within the Nitsche bounds. In case below the lower Nitsche bound we look 

at the remaining region A(r, p] = A \ in which = r=f. On the outer 

boundary of this ring we have h{x) = h°{x) = '"*|f|- On the other hand, it 
follows from the equation (13.5) that 

(13.29) l^^l = ^ r^{\h\) = ^v{r*) = 

\x\ \x\ 

Thus h is constant along each ray t , for r < t ^ p. This means that 
h{x) = r^^ = h°{x) in A\Q. 



CHAPTER 14 



Above the upper Nitsche bound, n ^ 4 



This chapter is devoted to a thorough discussion of the minimization 
problem when, in dimensions n ^ 4, the target annulus is conformaUy too 
fat. We shall see that the extremals are not radially symmetric. But first 
we need to examine related extremal problems for mappings on spheres. 



14.1. Extremal deformations of the sphere 

The study of extremal deformations of ring domains by using spheri- 
cal coordinates leads us to a variational problem for mappings of the unit 
sphere. This problem, to be explored later, is the following. Among all 
homeomorphisms $ : S"-^ §""1 of Sobolev class ^i'"(S"-i, S^-^) find 
the one which minimizes the energy integral 

(14.1) r[$] = / [a' + (n - l)p$]'] ^ 

where a is any given number. Here stands for the normalized Hilbert- 
Schmidt norm of the tangent map 

(14.2) : S"-i ^ Ty S'*-^ , y = $(x) 
That is, 

(14.3) [D^]^ = Tr [D*^D^] 

n — 1 

which equals 1 for the identity map. An obvious question to ask is whether 
the identity map 

(14.4) id : S"-^ ^ S"-^ 

is the minimizer. Naturally it is tempting to apply Jensen's inequality. 



(14.5) r[$] ^ 



a^ + {n-l)-f [D<^f 



In dimensions n = 2 and n = 3, one may appeal to Holder's and Hadamard's 
inequalities to find that 

2 2 

(14.6) / ^ (-[ [D^^-A ^ (-[ J{x, $) dx] = 1 
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14. ABOVE THE UPPER NITSCHE BOUND, n ^ 4 



As usual J(x, denotes the Jacobian determinant of the pullback via 
$ of the standard (n — l)-form uj on S"~^. 

(14.7) 4 J{x, $) dx = deg $ = 1 

Unfortunately, in dimensions greater than 3, the infimum of the integrals 



/ 



is not attained for the identity map. This infimum equals zero. Here is a 
computation which, in addition to illustrating this fact, provides a method 
for constructing more sophisticated examples. Consider a permissible map 
. §n-i _^ §"-1^ stretching a spherical cap of radius e around the north 
pole onto the entire sphere S""^.^ The rest of the sphere is shrank into the 
south pole. 



e _ 


S 




















N 



Figure 1. A spherical cap to be stretched around the sphere. 



(14.8) 



Elementary geometric considerations show that 

I otherwise 
whereas |S^| ^ e^"^. Hence, for every l^p<n— Iwe find that 

(14.9) / [D^'{x)Ydx = O{e''-^-P)^0 as e ^ 

the infimum being equal to zero. The above computation suggests that we 
have to express the integrand of T[$] , 

[a^ + (n- l)[i:»$]2]^ 
as a convex function in [D$]"^^. For this purpose we introduce 

' -1 H 

(14.10) F{s) = + (n- l)s^ ' , ^ s < oo 



is a weak W^'"-\\\mt of homeomorphisms of §" ^ onto itself. 
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Of particular interest to us will be the lower bound of F by a convex function 
F* such that 



(14.11) 



F* = F*is) ^ F{s) 

F*{1) = F{1) = (a^ + n- l)t 



14.2. Random variable setting 

It is both illuminating and rewarding to consider even more general 
setting of the variational integrals such as (14.1). We just replace 
by a general measurable function. 

Let (S,/x) be a probability measure space. We shall consider random 
variables X : S ^ M+ whose integral mean is at least one; 

(14.12) / Xdn^l 

Js 

The energy of X is defined by the formula 



(14.13) £[X] 



2 



+ (n- l)X"-i 



2 



where < a < oo and n = 2, 3, ... We look for the parameters a for which 
the constant function X = 1 is a minimizer. That is, 

n 

(14.14) M£[X]=£[1] = [a^ + n-l]^ 

This is certainly true for all a if n = 2, 3, because F is convex. However, in 
higher dimensions F changes concavity. It has an inflection point at 

n— 1 

/ In — R \ 

(14.15) So 

Precisely, F is concave for ^ s ^ So and convex for s ^ So- We shall first 
examine the case < a ^ a^, where the upper bound q:„ is determined by 
the equation 




(14.16) K + n-1) ^ {al-l)=a: 

We have 



(14.17) 1< a„ < 

Case 1. < a ^ a„. This means that 



n 



n 



n — 2 

(14.18) (a^ + n-l) ^ [a"^ - l) <^ a- 

It is important to observe that in this case the tangent line oi F = F{s) at 
the point s = 1 lies entirely below its graph, see Figure 2. 
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We then consider a convex lower bound F*{X) ^ F{X) which coincides 
with that of F for s ^ 1 and extends along the tangent line for ^ s ^ 1. 
Jensen's inequality yields 

S[X] = J F{X) dn^ J F*{X) dn ^F* (^J X dfi^ 

(14.19) ^ F*{1) = F(l) = S[l] 

Furhermore, equality occurs if and only if X = 1. 

Case 2. a > an- Thus the tangent line to F at s = 1 intersect the graph 
o{ F = F{s) at some point near the origin, provided a is closed to «„. This 
geometric observation suggests that we must look for the extremals which 
assume exactly two values. The best choice turns out to be when X assumes 

exactly two values, and j -To see this we split the sample space S 
into two parts 

(14.20) 

where 

(14.21) 

and define 



S = Si U S2 



m(Si 



n-l 



and iJ, (S2) = 1 



On 

a 



n-l 



(14.22) 



Xo 



{an) 



n-l 







on Si 
on S2 



The energy of Xo is, therefore, easily computed as 

„2- 



(14.23) 



£[Xo] 



+ {n- 1 



a 



at. 



/x(Si) + a"M(S2) 



Taking into account the definition of q;„ at (14.16), we arrive at the formula 
(14.24) f[Xo]=a" + 6a 
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where 



n (ai + n — l) 

(14.25) 6 = — !^ " ' 



n-2 
2 



There remains the question as to whether Xo possesses the minimum energy 
among all random variables of mean at least one. To this end we appeal to 
Lemma 9.1 in its borderline case when a = a(a„) = 1. By formula (9.4) we 
obtain 



(14.26) 



+ in- l)X"-i 



' ^ a'' + haX 



for every random variable X : S — > IR_|_. Equality holds if and only if X 

assumes exactly two values and (^^^ • Upon integration over S we 

conclude with the desired estimate 

(14.27) £[Xo\^a'' + ha 

For the equality to hold the distribution function of the random variable X 
must coincide with that of Xo . In particular, the constant function X = 1 
is not the extremal one. 

14.3. Pulling back a homothety via stereographic projection 



Before we see what can happen if a > y , let us consider an example 
of a non-isometry $ : S"~^ — >■ with 



/ 



in— 1 

§n-l 

Here is a construction of such a map. 

Let n : §"~^ R"^^ denote the stereographic projection of the unit 
sphere S"~^ C M" through the south pole onto R"~^. Given any positive 
number A we consider the homothety f = ■ defined by 

/'^(x) = \x for X G R""^. Conjugate to is a conformal mapping $ = 

$A . gn-1 _^ gn-1 

(14.28) = n-^ o o n 

We call ^ = the spherical homothety. Let x = (cos0,ssin^) G §"~^ be a 
point of longitude 5 G S""^ C M"~^ and meridian ^ ^ ^ tt, see Section 4.3. 
The south pole corresponds to ^ = tt. Thus $ : §"~^ — )• S"'~^ is longitude 
preserving. Let (p = (p{6) denote the meridian of ^{x) = {cos (p,s sin ip), see 
Figure 3. 

Geometric considerations give the following explicit formula 

(14.29) tan| = Atan^ 
see Figure 4. Hence 

(14.30) (p{e) = 2 lan~^ tan ^ 
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Figure 3. Spherical homothety $ : ^ ^ ^ via stere- 
ographic projection. 




Figure 4. Computation of the meridian of x and 



and 

^''■''^ ^(^^^ icos^| + Asin^i >" 

(^(0) = 0, (^(7r)=7r, (^(0) = A, (/j(7r) = ^ 

see Figure 5. 

Further differentiation shows that 

(14.32) m = 1(^1- >^ymf sine 

Thus (f is convex if < A < 1 and concave if A > 1. The extreme values of 
ip are assumed at the end-points, 

(14.33) min|A, ^1 ^ vp(6') ^ max|A,^ 
Implicit differentiation of (14.29) yields 
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A = 1 




TT B 



Figure 5. The meridian functions. 



Hence 
(14.35) 



sm if 
sin^ 



for every parameter A > 0. 

Next, we return to the mapping $ : S'*"^ — >■ S""-*^ and its tangent bundle 
map : TS"-^ ^ TS"-\ 



(14.36) 

Thus 
(14.37) 

In particular, 
(14.38) 



2 def 1 



n-1 



+ („ - 2) 



sin^ ip 
sin^e 



sin^ (f 



. . sm (p 
sm ^ 



mm 



Finally, let us entertain the reader with the following computation 
(14.39) = detD$^ = 

Hence we find that 



sm V \ " ^ sin" ^ 
sin^ 



sin"-2 Q 



4 

(14.40) 

as expected. 



. sm" {p 



.n-2 , 



sm 



sm" 
/J'sin'^-^^c/^ 



sm 



n-2 , 



de 
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14.4. Back to the variational integral T[^] 
Recall that 

(14.41) r[$] = / I + (ra - l)p$]2 |i 

for homeomorphisms $ : S"""^ — > S""-*^. We will test it with $ = so that 

1 



(14.42) 



n — 1 

Now we recall the function 



02 + (n-2)^ 



sin 



sin^ if 



(14.43) 



F{s)= a^ + (n-l)s"-i 



, s > 



The first and second derivatives of F are given by the formulas 



(14.44) F{s) = n + (n - l)s 
and 

(14.45) F{s) = n + (ra - l)s^ 



2 

n-1 



s"-i > 



ri — 3 



-a 



Hence F is convex if n = 2, 3. For n ^ 4 the function F is concave in the 
interval 



(14.46) 



^ s ^ 



n-l 

n-3\~ 



n — 1 



a 



n-l 



From now on, we consider the case 



n — 1 
n — 3 ' 



(14.47) a" > 
Given any parameter 

(14.48) 1< A ^ 



n — 3 
n-l 



We examine the spherical homothety $ = $^ : ^ — > S'* ^. First notice 
that 



(14.49) 



[D^ 



n-l 



n — 3 
n — 1 



■ a 



n-l 
2 



point-wise everywhere. Since [D^] is not identically equal to one, by con- 
cavity argument we conclude with the following strict inequality 



/ [a^ + {n - l)[D^^f = / 

J§ri-1 J§ri-1 
(14.50) < 



a 



+ (n-l) (p$]"-i)~ 



2 + (n - 1) (-[ [D<^r-A 
[a^ + n - 1] ^ 
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jn— 1 , sn— 1 



Thus the identity id : §""^ is not a mmimizer. 



14.5. The failure of radial symmetry, Proof of Theorem 1.11 

Throughout this section we make the following standing assumption on 
the moduli of A and A* 

(14.51) ModA^ModA* 

Theorems 1.9 and 1.14 tell us that in dimensions n = 2,3 the radial mappings 
are unique (up to a conformal automorphism of A) minimizers of both energy 



(14.52) <4 = \\Dhr, he^{A,A*) 

Ja 
f II Dh i " 

(14.53) = y^^^' hen^,A*) 

This is also true for n ^ 4 provided the target annulus A* is not too fat, as 
in Theorem 1.10 and Theorem 1.15. Here we construct examples to show 
that these results do not extend to the full range of moduli at (14.51). 

Theorem 14.1. Suppose 



/n — 1 

(14.54) ModA* > W -Mod A, n>4 

\ n — 6 

Then 

(14.55) mf / " , , , " > mf / " , , , " 

where ^(A , A*) stands for the class of radial stretchings h : A A*. 

Proof. In Section 8.2 we have shown that the infimum in the left hand 

side is attained for 

ModA* 



(14.56) h = h'^(x) = nr-''\x\'^-'x, a 
^ ' V / * I I > ModA 

The =^-energy of /i" is easily computed 

n 

(14.57) =^/,a = (a^ + n - 1) 2 ModA 

We now test the infimum on the right hand side of (14.55) with the spherical 
mapping (which is not radial), 



(14.58) /i(x) = r,r-°|xr-^$^ 



where <1>^ : ^ S" ^ is the spherical homothety defined at (14.28). If A 
is sufficiently close to 1, but difi"crcnt from 1, then 

f r -1 — 

(14.59) 4 + {n-l)[D^^f ^ <{a^ + n-lY 

J§n-1 L J 
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see (14.50). Now, the ^-energy of h is computed as follows 



\hN\' 



+ (n-l 



\M 
\h\' 



2\ 2" 



= / 

J A 



M2\ 2 



dx 



= 4 fa2 + (n-l)[D$^]2)' -Mod. 

J§n-1 \ J 

71 

< (a^ + n- 1) 2 ModA = =^fta 



(14.60) 
as desired. 

The spherical homothety : S"~^ — >■ S"~^ works as well for the n- 
harmonic energy, though a computation is more involved. It results in the 
proof of Theorem 1.11. 



Theorem 14.2. Let 4 an. 
' \Jn — 1 + a/tT 



(14.61) 8n = 



_ 1 
3\ 2 



exp 



n 



n\/n — 1 



tan ^ — 3 



^ ^/n. 



Consider the annuli A = A(r, R) and A* = A(r=K, -R^,), such that 
(14.62) 
Then 



^ ^ ^ g-('^n) 

l<7<<^n and ->J^^^ 



(14.63) 



inf / 

fee^(A, A*) 



> inf 



mil" 



Proof. First we show that the infimum on the left hand side is attained 
at the radial n-harmonic map 

X 

\x\ 



(14.64) 



h{x) = H{\x\)-- 



of the form H{t) = \H^(kt), for suitable parameters k > ^ and A > 0. 
Indeed, the energy of radial mappings is given by 

rR r 



(14.65) 



f 

J r 



t^H^ + (n - 1)H^ 



i dt 
t 



By standard convexity arguments the infimum is attained and the minimizer 
satisfies the Lagrange-Euler equation. This equation simply means that h is 
n-harmonic. Since Mod A ^ Mod A* we are in a situation in which h lies in 
the class generated by the principal solution H^; that is, H{t) = \H-{kt), 
as claimed. These parameters are uniquely determined by the size of the 
annuli A and A*. Precisely, we have 



(14.66) 



XH^{kr) = r* and XH^{kR) = i?* 
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We eliminate the factor A by dividing the above equations, and then find k 
from the equation 

To find k we observe that the function Q{t) = , ^ < t < cjo, is 

increasing from to . Indeed 

tQ{t) trH(tr) tRH(tR) , , . „x . 

because the function ij- = r/_(s) is decreasing, see Figure 8 in Chapter 5. 
This latter statement is immediate from equation (5.29); that is, 

1 + (ryg.-!) = ) where ii"(s) = ii"_(s) , for simplicity 

Here — c > and |i?(s)| is increasing. Also note that 

tR 

see (5.59). Now the equation (14.67) has exactly one solution. Our hypoth- 
esis at (14.62) is equivalent to 



(14.68) g(fco) = ^<Q{^ 

-Tt* \ -ft 

Since Q is increasing we infer that koR < 6n- Next, since rf- is decreasing 
then foT r < t < R, we have 

rj^{kot) > ri^{koR) ^ ri-{dn)- 
Now we appeal to formula (14.61) and (5.47). Accordingly, 

'n — 3' 

On the other hand, in view of (5.53) and (5.49) 



1 1 In - I 



ui6n) r-i(5„) Vn-3' 
Therefore, r]-{kot) ^ for every r < t < R. This means that 



(14.69) r/^ (t) > J - — \ , for all r ^ i ^ i? 
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Now, we return to the computation of the infimum in the left hand side of 
(14.63) to obtain 



inf 



Dh\ 



R 



t^H"^ + (n - l)H' 



f dt 



(14.70) 



'■R n rif 



Then we test the infimum in the right hand side of (14.63) with the mapping 



(14.71) 



hxix) = H{\x\) 



where, as in the previous case, : S" ^ — >■ S" -"^ is the spherical homothety. 
We find that 



inf [ \\Dh\r ^ [ WDhxir 

^(A,A*)A J A 



R 



mr-f \viit) + {n - i)\D^'\' 



f dt 



(14.72) 



dt 



< un-i / [H{tr[r^l{t) + n-l]^j 



inf / WDhr 



The above strict inequality follows from (14.50) applied to a = ri^(t) > 



where we have chosen A / 1 sufficiently close to 1. 



CHAPTER 15 



Quasiconformal mappings between annuli 



In this final chapter we present an apphcation of free Lagrangians to 
obtain sharp estimates for quasiconformal homeomorphisms /i : A — )■ A* 
between annuli A = A{r,R) and A* = A(r*,i?*). With the aid of Mobius 
transformations (reflections about the spheres and (n — l)-dimensional hy- 
perplanes) we may assume that h preserves the orientation and the order 
of the boundary components. We shall employ the operator norm of the 
differential matrix, commonly used in the literature [29, 43, 44]. Accord- 
ingly, a homeomorphism h : A A* of Sobolev class ^^^^''"(A, A*) is 
quasiconformal if 

(15.1) \Dh{x)\'^ ^KJ{x,h) a.e. in A 

The smallest such number K ^ 1, denoted by Kq = Koih), is called the 
outer dilatation of h. The inner dilatation is defined to be the smallest 
number Kj = Kj{h) ^ 1 such that 

(15.2) |D«/i(x)|" ^ J(x,/i)"-^ a.e. in A 
Note the relations 

(15.3) Ko ^ K^~^ and Kj ^ K^,'^ 
For the power stretching h{x) = \x\°'~^x, a > we have 

(15.4) Ko = max {a~^ , a""^} and Kj = max {a^"" , a] 

which shows that both estimates at (15.3) are sharp. Since J{x, h) is inte- 
grable and X < oo we see that h G >r^''*(A, A*). 

Theorem 15.1. Suppose h : A —> A* is a quasiconformal map between 
annuli. Then 

1 /ModA*\""^ 

This estimate is classic in the theory of quasiconformal mappings, see 
the pioneering work by F.W. Gehring [16]. In the proof below wc shall not 
appeal to any advances in Quasiconformal Theory or PDEs. In fact, our 
proof provides a method to tackle the uniqueness problem in the borderline 
cases of (15.5), which seems to be unknown in higher dimensions. 
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Theorem 15.2. If one of the two estimates at (15.5) becomes equality, 

then it is attained only on the corresponding extremal mappings h : A 
A* of the form 



(15.6) h{x) = ^/r^ ( -j^ I ^' ( — ) ' r ^ ^ ii: 

where 




(15.7) a 



Mod A* 
Mod A 



K^-' if Mod A* ^ Mod A 



L K/-" if Mod A* ^ Mod A 



Here the spherical part $ : ^ ^ can be any homeomorphism 

satisfying 

(i) volume condition: J{u}, #) = ±1 , for a.e. to E E>'^~^ 

(ii) a- contraction condition: meaning that 

• I D^{uj) I ^ a , for the equality in the upper bound of (15.5) 

• I {D<^{oj))~^ I ^ 1/a , for the equality in the lower bound of (15.5) 

Note, that the only volume preserving homeomorphisms $ : §^ 
are isometries. Thus, in dimension ri = 2, the extremal quasiconformal 
mappings take the form 

(15.8) h{z) = \\z\'^-^ z 



where A is a complex number of modulus r r^ Example 15.4 shows that in 
higher dimensions there exist non-isometric homeomorphisms <I> : S"""^ 
§"~i which are volume preserving and satisfy the a-contraction condition, 
see also [42]. The exception is the conformal case of a = 1 for which the 
contraction condition |Z)$(a;)| ^ 1 together with J(a;, $) = ±1 imply that 
$ is an isomctry. As a corollary to these observations we obtain Schottky's 
theorem (1877), [45], in M". 

Theorem 15.3. An annulus A = A(r, R) can be mapped conformally 
onto A* = A(r*,i?*), if and only if y = Moreover, modulo isometry 
and reseating, every conformal mapping takes the form 

(15.9) h{x) = 



the identity 
the inversion 



For both Theorems it involves no loss of generality in assuming that h 
preserves orientation and the order of boundary components of the annuli. 
And we do so from now on. 
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15.0.1. Proof of Theorems 15.1. Let us first prove the inequahty at 
the right hand side of (15.5). We begin with the identity (7.6) and, after 
Holder's inequality, use the free Lagrangian (7.4) 



d\h\ A *d|x| f /r^,n* h X \ dx 



1^1 kl""-^ Ja\ \h\'\x\/ |/i| |a:|"-i 

1 n-l 

\Dh\dx ( f " / f dx 



Mod A* = 



^ Ja \h\ |x|"-i ^ VA l/il" J KJa \x 
= Kg (Mod A*)- (ModA)"^ 

as desired. 

For the left hand side at (15.5) we begin with the identity (7.12) to 
compute in the similar fashion that 

J A kl JaY VI \h\/ \x\\h\^-^ 

r \D^h\dx ( r \Dih\^\~ f f dx 

JaW¥^^[Ja 



\h\^ / \Ja \x 



(15.11) ^ Kf (^jTZ^g^^) - (ModA)^ 



\h\- 

1 n-l 

(ModA*)~ (ModA)' 



as desired. 



15.0.2. Proof of Theorem 15.2. Concerning the borderline cases, let 
h : A A* denote the extremal mapping for the right hand side of (15.5). 
We simplify the matters by assuming that r = r* = 1. Note that the power 
mapping 

ModA* 



(15.12) /i"(x) = br-^x, a 
^ ^ V y I I > ModA 

is among the extremals. The borderline equation 

(15.13) (ModA*)'*-^ = Ko (ModA)"-^ 

implies that we have equalities everywhere in (15.10). This amounts to the 
following two conditions. 



0)1^*^-M'r)= = ' « - a positive 

constant 

and 



(ii) \Dh\'' = KoJ{x,h) 



(15-18) D*h-j-^ = a- 
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The first set of equations is fulfilled if and only if 

h X \h\ X 

(15.14) D*h ■ — = \Dh(x) \ — = a\-j — , a - a positive constant 

\h\ \x\ \x\ \x\ 

We identify a from the following equations 

(15.15) ^J^^Ko'-^ 

Upon integration over the annulus A we obtain 

(15.16) a" Mod A = Mod A* 
which in view of (15.13) yields 

Now, an interesting nonlinear PDE arises from (i) and, in a more direct 
manner, from (15.14) 

h _ X 

Rather unexpectedly we can easily solve this equation for \h\. Let us express 

it as 

(15.19) Vlog|/i|^ = QVlog|a;|2 

Hence, in view of the normalization \h{x)\ = 1 for |a;| = 1, we find that 

(15.20) \h{x)\ = 
Now, h takes the form 

h{x) = |a;|"^'(x) , with ^' : A ^ S""^ 

The map h, being a homeomorphism in /^^'"(A, A*), is differentiable almost 
everywhere, then so is ^. From (15.14) we also find the operator norm of 
the differential 

(15.21) \Dh{x)\ = = «kr"^ = \Dh"ix)\ 
This, in view of (ii), yields 

(15.22) J{x, h) = a|a;|"°^-" = J(x, /i") 

Further examination of the mapping ^ : A — > S""^ will reveal that its 
normal derivative vanishes almost everywhere. To this end, let a;o € A be a 
point of differentiability of ^. For small real numbers e we consider Taylor's 
expansions of order one, 

(15.23) ^{xo + exo) = ^(xo) + e*iv(a;o) + o(e) 
and 

(15.24) h{xo + eXo) — h{xo) = e[Dh{xo)]xo + o(e) 
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Upon elementary computation, letting e go to zero, we arrive at the following 
estimate 



(15.25) \a^{xo) + *iv(xo)| ^ a 
We square it to obtain 

(15.26) + 2a(*, ^-jv) + l^'jvl^ ^ 



It is important to observe that \I' is orthogonal to its directional derivatives, 
because = |*|^ = 1. Thus, in particular, "^^(xo) = 0. This simply 

means that 



(15.27) *(x) = $ ( ^ ) , where $ : S'*"^ ^ ^""'^ 

Obviously, $ e /r^'"(S"-\ S""^) is a homeomorphism of S""^ onto itself. 
It induces the linear tangent map 



(15.28) D^iio) : T^S""^ ^ r^^S""^ , a = $(a;) 

The Jacobian determinant of with respect to the standard volume form 
on will be denoted by 



(15.29) J(a;,$) = detD$ 

It relates to the Jacobian determinant of h by the rule 

(15.30) J{x,h) = a|x|""-"($,$r2 X ... X $T„) 

(15.31) = a|xr-"J(a;,$), ^ = ^ 

Thus, in view of (15.22) we see that $ : §"~-^ — >■ §"~^ is volume preserving. 

Here we can take the cross product of directional derivatives with respect 
to an arbitrary positively oriented orthonormal frame T2, ...,r„ in T^^S""^. 
Thus, in view of (15.22), we obtain 



(15.32) J(a;,$) = |$t2 x ... x $t„| ^ I^Tal • • • |^t„| ^ 



\n-l 
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Next, we look at the Cauchy-Green tensor of h in terms of the orthonormal 
frame AT, T2, r„ at T^A 



D*hDh 



- - Hn - - 

- - - - 

- - hT„ - - 



hN ■ ■ ■ hT„ 



(15.33) 



|2o-2 



- — — 

- $T2 - - 

- *T„ - - 



a$ $T2 • • • ^T„ 



= \x 



\2a-2 



a' 








and 
(15.34) 



Let ^ A2(a;) ^ ... ^ An(a;) denote the singular values of D^, meaning that 
A|(x) ^ ... A^(a;) are eigenvalues of D*^D^. It follows from (15.10) that 
the numbers 

A2(x)|x|°-\...,An(ar)|x|°-^ and a|a;|"-^ 

are the singular values of Dh. Since q;|x|"^~''^ = \Dh\, this latter number is 
the largest singular value, 

(15.35) ^ A2(a::) ^ ... ^ A„(x) ^ a 
In particular, 

(15.36) [D<P{x)] = Xn{x) ^ a 
as desired. 

This computation also shows that the inequality |D<I'(x)| ^ a is both 
sufficient and necessary for the equation \Dh{x)\ = a|x|"~^. 

A backwards inspection of the above arguments reveals that h{x) = 

(y'^l) extremal map if and only if |I?$(x)| ^ a and J{u, $) = 1. 

It is rewarding to look at the inverse map f = : A* A, which 
takes the form 



(15.37) 



m = \y\'-^-' (^) 
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We observe that for every quasiconformal map h and its inverse h^^ = f it 
holds 

(15.38) Ko{x, h) = Ki{y, /) , y = h{x) 

Now in much the same way we find all extremals for the left hand side of 
(15.5). Theorem 15.2 follows. 

Let us finish this subsection with an example of the volume preserving 
mappings of C M^. 

Example 15.4. For every e > there exists a homeomorphism $ G 
>r^'°°(S2,S2)^ not rotation, such that - uj\ ^ e and J(a;,$) = 1 for 

almost every w G S^. 

15.0.3. Construction of We shall view $ as a small perturbation 
of id : ^ S^. The construction of such perturbation is made in three 
steps. In the first step we project the sphere + + = 1 with two 
poles removed, z ^ ±1 onto the cylinder a;^ + y^ = l, — l<z<l. The 
horizontal rays from the axis of the cylinder project the points of the sphere 
towards the surface of the cylinder. This projection, sometimes attributed 
to Archimedes, is well known as Lambert's Cylindrical Projection (Johann 
H. Lambert 1772). In the second step we cut the cylinder along a path 
from the top to the bottom circles and unroll it flat. The third map is a 
piece- wise linear perturbation of the triangle ABC inside the fiat region, see 
Figure 1. It keeps A, B, C fixed while permuting X, Y, Z in the cyclic way, 
X— T-Y— T-X. It is geometrically clear that such piece- wise linear 
deformation preserves the area. Moreover, choosing X, Y, Z close to the 
barycenter makes the deformation arbitrarily close to the identity. Finally 
we roll it back onto the cylinder and use inverse of Lambert's projection to 
end up with the desired area preserving perturbation of the identity. 

Remark 15.5. As pointed out by the referee, a simpler example can be 
given : (9,if) i— )■ (9 + g(9),(p), where 6 ,if are the longitude and latitude, 
and g is any Lipschitz function. 

Acknowledgements. We thank Leonid Kovalev for valuable discus- 
sions on the subject of this paper and the referee for very careful reading of 
the paper and many useful suggestions. 
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Figure 1. Constructing an area-preserving Lipschitz homeomorphism. 
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